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ABSTRACT 


The  present  work  investigates  the  aerothermal  effects  on  the 
electro-optical  performance  of  a  model  infrared  seeker  missile  at 
Mach  4  and  small  angle  of  attack.  Two  window  dome  configurations  are 
studied:  an  optical  hemispherical  and  ellipsoidal  shape  window 
attached  to  a  3:1  ogive  nose.  A  three-dimensional  thin-layer 
compressible  flow  solver  is  used  to  compute  the  flow  field  in  front 
of  the  window  and  the  window  aerodynamic  heating.  Grid  adaptation  is 
performed  in  order  to  accurately  compute  the  flow  field  and  capture 
the  bow  shock.  Noise-Equivalent-Temperature-Difference  degradation 
due  to  "background"  noise  of  the  hot  window  is  derived  and  computed 
for  3  levels  of  array  non-uniformity  for  a  range  of  mean  dome 
temperatures  covering  the  whole  supersonic  regime.  Acquisition  range 
for  a  staring  sensor  is  estimated  for  several  maritime  scenarios 
modeled  by  LOWTRAN  code.  Background  photon  flux  distribution  on  an 
array  generated  by  the  hot  dome  is  computed  using  radiation  transfer 
methods  and  the  computed  dome  temperature.  It  is  found  that  array 
non-uniformity  has  a  strong  influence  on  the  seeker  performance. 
This  study  indicates  that  the  non-unif ormity  of  the  dome  temperature 
field  has  a  significant  influence  on  the  array  fixed  pattern  noise 
induced  by  dome  emission  and,  as  a  result,  on  the  seeker  detection 
and  tracking  ability. 
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I.  INTRODUCTION 


Flight  at  high  supersonic  speeds  through  the  atmosphere  is 
characterized  by  high  aerokinetic  heating  rates  and  gas  compression 
in  front  of  the  missile.  Large  density  variation  and  heating  of  the 
window  significantly  affects  the  performance  of  a  vehicular  infrared 
(IR)  sensor,  leading  to  major  design  problems  of  IR  seeking 
missiles.  Presently,  IR  windows  that  operate  in  the  low  supersonic 
flight  regime  exist  and  are  in  operation.  However,  future  IR  and 
dual  mode  tactical  missiles  will  require  IR  systems  which  are 
capable  of  effective  operation  in  much  more  severe  flight 
environments.  The  durability  and  survivability  of  the  dome  plays  a 
major  role  in  defining  the  flight  envelopes  of  these  missiles. 

For  the  case  of  a  classical  IR  seeker  installed  at  the  tip  of 
the  missile  nose,  the  outermost  effect  of  the  flow  field  on  the 
optical  performance  is  due  to  the  large  density  jump  across  the  bow 
shock  wave  in  front  of  the  IR  window.  This  density  gradient  refracts 
optical  rays  and  acts  as  a  lens.  Gas  behind  the  shock  is  highly 
compressed  and  heated.  Density  variations  in  this  compressed  flow 
downstream  of  the  shock  create  refractive  index  variations  and  have 
the  effect  of  a  gradient  index  lens.  The  principal  optical  effects 
are  resolution  (blur  and  defocus)  and  boresight  error.  Another 
effect  of  the  compressed  flow  is  aerodynamic  heating  of  the  window 
itself.  As  the  IR  dome  heats  up,  its  self-emission  may  adversely 
impact  the  operation  of  the  seeker  system  in  a  high-speed  flight 
situation.  The  added  background- induced  noise  will  lower  the 
sffsctive  performance  of  the  system  in  terms  of  detection  range. 
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resolution  and  tracking  capability.  High  resolution  is  desirable 
for  detection  (including  clutter  rejection),  tracking  accuracy,  and 
recognition.  The  boresight  error  can  affect  performance  in  several 
ways.  Boresight  fluctuations  (caused  by  turbulence,  shock 
interaction,  or  coolant  mixing)  can  blur  images  and  induce  angle 
noise  in  the  tracking  system.  Changes  in  boresight  caused  by  body 
motion  (e.g.  changing  shock  position  or  window  distortion)  appear  as 
position  errors  or  tracking  rate  errors  that  adversely  affect 
sensors  used  in  guidance  and  navigation. 

Significant  radiation  from  the  compressed  gas  comes  from  heating 
by  its  passage  through  a  strong  shock  wave.  This  effect,  however,  is 
only  significant  at  higher  Mach  numbers  M  >1) ,  when  the  perfect 
gas  assumption  breaks  down  and  chemical  reactions  must  be 
considered. 

A  mission  analysis  [1]  for  the  future  tactical  missile 
generation  revealed  a  requirement  for  high  supersonic  speeds  of 
operation.  This  will  produce  very  high  aerodynamic  pressures  and 
thermal  stresses  on  the  optical  window  that  protects  the  IR 
detection  system.  Severe  environmental  conditions  encountered  at 
very  high  aerodynamic  speeds  has  forced  the  development  of  new 
exotic  materials,  capable  of  withstanding  such  high  temperatures  and 
erosion.  Traditionally,  the  shape  of  most  IR  windows  or  domes  that 
protect  the  IR  seeker  has  been  a  truncated  hemisphere.  Although  this 
configuration  has  some  advantages  from  the  optics  point  of  view, 
because  of  easier  fabrication  and  design,  it  imposes  a  considerable 
aerodynamic  drag  penalty  on  the  vehicle  [2].  Therefore,  an 
ellipsoidal  window  is  proposed.  Such  window  design  would  decrease 
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the  aerodynamic  drag  of  the  vehicle,  increasing  its  range.  In 
addition,  it  is  expected  that  it  will  delay  boundary  layer 
transition  and  reduce  aero thermodynamically  imposed  stresses  within 
the  IR  window  material  [3] .  As  a  result,  the  self-emission  photon 
radiation  from  the  window  into  the  detector  will  decrease  [2]  [4], 
and  the  detection  range  will  improve. 

The  temperature  field  on  the  window  is  obtained  nvimerically, 
using  the  OVERFLOW  flow  solver.  Comparison  with  available  data  was 
conducted  to  validate  to  some  degree  the  computation,  and  to  give 
confidence  in  the  temperature  results  obtained.  Using  this  data,  the 
aero-thermal  effects  which  affect  the  imaging  performance  of  an  IR 
seeker  were  evaluated. 


RESISTS  RAIN 
EROSION 


Figure  1-1.  Infrared  Dome  Design  Constraints 
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II.  BACKGROUND 


A.  BLUNT  BODY  STRUCTURE 

A  typical  blunt -body,  supersonic  flow  field  is  shown  in  Fig. 2-1. 
This  figure  shows  the  general  structure  of  the  flow  field  in  front 
of  the  window  mounted  at  the  nose  of  a  supersonic  vehicle.  The  shock 
wave  is  detached  and  normal  shock  relations  may  be  used  to  estimate 
the  change  in  density  along  the  stagnation  streamline.  A  complete 
flow  field  and  shock- structure  solution  depends  on  body  geometry  and 
flight  Mach  niomber.  Such  solutions  are  complex  because  of  mixed 
subsonic  (M<1)  and  supersonic  {M>1)  regions  in  the  flowfield.  The 
governing  equations  are  mathematically  elliptical  in  the  subsonic 
region  and  hyperbolic  in  the  supersonic  region.  Time -dependent 
numerical  techniques  are  used  to  determine  flow  fields  surrounding 
blunt  bodies [5].  An  empirical  formula  for  the  shock  stand-off 
distance.  A,  at  the  tip  of  a  conical  body  with  a  hemispherical  nose 
is 


A/R  -  0.143e 


[3.24/M^] 


(2-1) 


where  R  is  the  nose  radius .  This  equation  shows  that  the  shock 
standoff  distance  decreases  with  Mach  number.  In  this  formula  it  is 
assiomed  that  the  detached  shock  wave  is  a  hyperbola  that  is 
asymptotic  to  the  free  stream  Mach  angle .  The  equation  for  the 
coordinates  of  the  shock  around  a  hemisphere-cone  [6]  is 


X  .  R.A-R  cotan^Q 

C 

\ 

^  y^tan^Q 

1/2 

-1 

(2-2) 
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Figure  2-1.  Illustration  of  the  flow  field  in  front  of  the  missile's  nose. 

B.  AERODYNAMIC  HEATING  AND  HEAT  TRANSFER 
1 •  The  Heat  Transfer  Process 

A  vehicle  operating  at  supersonic  conditions  is  heated  by  heat 
generated  within  the  boundary  layer  that  surrounds  it.  Heat  is 
produced  by  the  conversion  of  kinetic  energy  through  flow 
deceleration  in  the  boundary  layer,  as  shown  in  Fig. 2-2.  Heat  is 
transferred  to  the  body  surface  by  conduction  which  is  represented 
by  Fourier's  law  of  heat  conduction  given  by  [7] 
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q  =  -k 


(2-3) 


where  the  temperature  gradient  in  the  boundary  layer  is  evaluated  at 
the  wall . 


Figure  2-2.  Thermal  boundary  layer  temperature  distribution 
near  insulated  wall  [8J. 

Thermal  diffusion  by  conduction  is  a  complex  three-dimensional 
problem.  The  solution  of  this  problem  requires  the  formulation  of  an 
energy  balance  which  considers  the  thermal  energy  stored  at  each 
point  in  the  structure  as  well  the  thermal  energy  conducted.  In 
rectilinear  coordinates  this  energy  balance  is  given  as  follows: 


dT 

dt 


ot 


d 


dx  ‘ 


d 

dy^ 


d 


(2-4) 
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Figure  2-3.  Effect  of  heating  and  cooling  on 
temperature  distribution  [8J. 


Figure  2-4.  Distribution  of  temperature  near 
insulated  wall  [8J. 


In  practice  heat  transfer  is  computed  using  finite  differences  and 
finite  element  methods  for  heat  conduction.  At  elevated  temperatures 
a  surface  will  lose  heat  to  its  surroundings  by  radiation  according 
to  the  Stefan-Boltzmann  law,  given  by 


-  sarr;  -  (2-5) 

where  cp=5 . 67x10  ®  Wm  ^  K  ^  and  s  is  the  surface  emissivity  having 
values  between  0  and  1  depending  on  the  material,  is  the 
background  temperature  at  which  the  background  radiates  back  to  the 
surface.  There  is  also  the  possibility  of  radiating  from  the  skin  of 
the  vehicle  to  its  own  internal  structure.  While  the  form  of  the 
Stefan-Boltzmann  law  is  the  same,  the  temperatures  of  the  radiating 
surfaces  and  non-gray  absorption  and  emission  must  be  considered. 
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EXTERIOR 


2 .  Recovery  Temperature 

It  is  important  to  note  that  in  gases  the  adiabatic  wall 
temperature  ,  is  always  less  than  the  free-stream  stagnation 
temperature,  T^.  This  condition  is  in  contrast  to  what  one  may 
expect  from  the  well-known  temperature-Mach  number  relation 

1.  [  (y-l)  /2]M^ 

’  1.  [  (y-1)  /2]M^  ' 


viscous  effects  within  the  boundary  layer  decrease  the  velocity 
from  the  freestream  value  at  the  outer  edge  of  the  boundary  layer  to 
zero  at  the  surface.  Pure  compression  heating  would  give  a 
stagnation  temperature  of 

I  \ 


T  .  T 

OO 


(2-7) 


at  the  surface  (that  is,  in  Eq.  2-6  let  Tj  =  and  =  T„.  Then 
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is  the  value  of  temperature  where  the 


M,  =  0  and  M2  =  MJ  . 
velocity  is  zero,  i.e.,  at  the  stagnation  point.  On  the  other  hand, 
in  continuum  flow  the  velocity  at  the  wall  is  zero,  but  as  already 
stated,  the  adiabatic  wall  temperature  is  always  less  than  the  free- 
stream  stagnation  temperature 


Mach  Number 


Figure  2-6.  Stagnation  temperature  versus  Mach  number  at 
several  altitudes [9]. 

In  practice,  heat  is  conducted  outwards  through  the  boundary 
layer.  Viscous  dissipation  produces  additional  heat  and  the  actual 
temperature  reached  at  the  surface  (assuming  no  heat  transfer  into 
the  surface)  is  called  the  recovery  te7nperature[8]  .  Therefore,  for 
practical  computations,  it  is  convenient  to  define  a  quantity  known 
as  the  recovery  factor,  r,  which  is  a  measure  of  the  fraction  of  the 
local  free-stream  dynamic -temperature  rise  recovered  at  the  wall: 


10 


It  is  obvious  that  if  the  entire  plate  is  adiabatic,  energy 
considerations  dictate  that  the  mean  stagnation  temperature  with 
respect  to  mass  flow  must  be  equal  to  the  free-stream  stagnation 
temperature  T^„.  Since  the  wall  temperature  is  less  than  T^,  it 
follows  that  the  distribution  of  stagnation  temperature  within  the 
boundary  layer  must  be  of  the  form  shown  in  Fig.  2-4,  with  some 
portion  of  the  boundary  layer  having  a  stagnation  temperature 
greater  than  T^. 

The  value  of  the  recovery  factor  depends  upon  Prandtl  number  Pr 
which  is  a  ratio  measuring  the  effects  of  conduction  and  dissipation 
mentioned  above.  The  Prandtl  nximber  is  defined  by 


Pr 


— ~  ~  0.71  for  air  at  N.P.T. 
k 


(2-11) 


where  /j  is  the  coefficient  of  viscosity  and  k  is  the  thermal 
conductivity.  These  coefficients  are  strongly  temperature  dependent 
and  graphs  showing  this  dependence  are  drawn  in  Fig. 2-7.  It  can  be 
shown  that  the  recovery  factor  is  approximated  in  terms  of  the 
Prandtl  number : 


r  .  Pr  "  (2-12) 

where  n=l/2  for  laminar  flow  and  n=l/3  for  turbulent  flow.  While 
this  is  a  useful  approximation,  the  recovery  is  also  a  function  of 
the  pressure  gradient  and  local  Mach  number. 

The  adiabatic  wall  temperature  is  then  given  by  an  equation  such 
as 
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T  . 

e 


2C 


(2-13) 


\  p/ 

where  is  the  velocity  at  the  edge  of  the  boundary  layer  and  Cp  is 
the  specific  heat. 

In  missile  flight  the  dome  is  subjected  to  aerodynamic  heating 
resulting  from  the  convective  heat-transfer  process  described  by  the 
equation 


■  ITF  (2-14) 

r  w 

where  is  the  wall  temperature  and  the  recovery  (adiabatic)  wall 
temperature.  The  adiabatic  wall  temperature  is  that  temperature  at 
which  the  heat-transfer  to  the  wall  is  zero,  i.e.,  where  the 
temperature  gradient  in  the  boundary  wall  is  zero,  as  shown  in 
Fig. 2-3.  The  heat-transfer  coefficient  is  the  ratio  of  the  heating 
rate  to  the  thermal  driving  potential . 

The  heat  transfer  coefficient  h  in  Eg. (2-14)  varies  considerably 
with  flight  conditions.  It  is  often  more  convenient  to  work  with  a 
non-dimensional  heat  transfer  coefficient  called  the  Stanton  niimber 
defined  by: 


St 


p  u  c 

e  e  p 


(2-15) 
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Stanton  number  and  skin  friction  coefficient  are  linked  by 
the  Reynolds  analogy 


St 


C£ 

2 


(2-16) 


where  Q  is  the  skin-friction  coefficient  and  S  is  the  analogy 
factor. 


3 .  Heat  Transfer  to  a  Hemisphere 

The  aerodynamic  flow  over  a  hemispherical  nose  is  laminar  at 
the  stagnation  point.  As  the  air  accelerates  in  passing  from  the 
stagnation  point  to  the  base  of  the  hemisphere,  it  can  go  through 
transition,  and  the  boundary  layer  generally  becomes  turbulent. 
This,  happens  at  Reynolds  numbers  around  4x10^  [2]. 

For  the  laminar  case,  Sibulkin  has  shown  that  locally  the 
Stanton  number  on  a  hemisphere  stagnation  point  can  be  written 


St  .  0.  763 


<  r,  ^ 

u 


pu^D 
\  ; 


(2-17) 


where  D  is  the  diameter  •  of  the  hemisphere,  0  is  the  velocity 
gradient  and  /j  is  the  viscosity. 

All  of  the  fluid  properties  are  evaluated  at  the  local 
conditions  just  outside  the  boundary  layer.  In  the  case  of 
stagnation  this  corresponds  to  stagnation  temperature  and  Pitot 
pressure.  This  reduces  Eq. (2-17)  to[2] 
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St  . 


9.37 


0.38 


(2-18) 


y/u/D 


0.25 

\  mA 

[  5) 

where  {Cpp/k^)  ^  is  the  stagnation  Prandtl  number  (approximately 
0.72  for  normal  flight  conditions)  ,  p^uD/jn^  is  the  free-stream 
Reynolds  niimber  based  on  hemisphere  diameter.  The  subscript  o  refers 
to  total  or  stagnation  values. 

From  Ref  [1] , 


St  -  0. 763 


1/2 

\ 

P^uD 

■1/2 

(  \ 
C  U 

-0.6 

1/2 

(  \ 

V  ^  J 

J 

k 

\  ”  / 

s.p. 

'’-J 

i 

(2-19) 


The  subscript  0  denotes  the  total  or  stagnation  value  after  the 
air  passes  through  the  shock  wave.  If  it  is  assumed  that  /j  and 
the  perfect  gas  equations  are  used  then  Eq.  (2-18)  reduces  to  the 
general  laminar  Stanton  number  equation.  This  gives  the  heat- 
transfer  rate  at  the  stagnation  point  where  the  recovery  factor,  r. 


h,,sf2R  -  1 . 7c^^f^JIXM^  (1-0. 2Mf)  ^  (2-20) 

equals  1.  The  laminar  heat-transfer  rate  decreases  away  from  the 
stagnation  point .  At  the  stagnation  point  another  semi-empirical 
formula  can  be  used  to  calculate  the  heat  transfer [4] 

From  this  formula  C .  Klein  showed  that  the  most  severe  thermal 
environment  in  terms  of  temperatures  does  not  correlate  with  the 
highest  thermal-flux  situation [4] . 

When  the  boundary  layer  is  turbulent,  the  Stanton  number  has  a 
maximum  at  about  40°  from  the  stagnation  point,  then  decreases 
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further  aft. 


As  shown  from  the  Eq. (2-7)  at  supersonic  speeds  the  temperatures 
reached  are  quite  high,  which  raises  the  issue  of  suitability  of 
certain  materials  as  IR  transmissive  materials  in  a  high  temperature 
environment.  However,  it  must  be  determined  first  whether  the  time 
frames  are  such  that  the  surface  temperatures  can  approach 
stagnation- type  limits.  For  that  purpose,  one  may  consider  the 
dimensionless  Biot  n-umber 


Bi  . 


(2-21) 


which  characterizes  the  transfer  of  heat  across  the  boundary  layer 
in  relation  to  the  dome's  ability  to  transport  heat.  Here  L 
represents  the  dome  thickness.  If  the  wall  is  thermodynamically 
thin,  i.e.,  Bi«l  and  the  radiation  and  conduction  losses  are 
negligible,  then  the  aerodynamic  heat-transfer  rate  can  be  equated 
to  the  rate  of  enthalpy  increase  within  the  wall,  or 


h  (T  -  T  )  -  phc 

where  p„  represents  the  window  density  and  Cp„  the  heat  capacity  of 
the  window. 

If  there  is  a  step  change  in  the  aerodynamic  conditions,  e.g. 
velocity  (thus  a  step  change  in  h  and  T^)  ,  the  solution  is 


dT 

W 

Ht 


(2-22) 
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T  -T 

r 

T  -T. 

r  1 


where  equals  the  initial  temperature  at-  t  =  0 .  At  higher  Mach 
numbers (higher  window  temperatures),  the  radiation  heat  transfer 
plays  a  significant  part.  If  this  is  to  be  considered,  the  heat- 
balance  equation  becomes 


h(T-TJ  -  pic 


.  oe(T)T„  . 


(2-24) 


The  Eq. (2-24)  can  be  solved  for  a  series  of  constant  temperatures  to 
yield  a  more  accurate  temperature  history.  The  true  equilibrium 
temperature  (  dT/dt=0  )  can  also  be  found  from  the  same  equation. 


Mach  Number 

Figure  2-8.  Time  constant  for  aerodynamic  heating [2]. 
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Eq. (2-25)  is  used  to  give  an  indication  of  the  relative  rate  of 


temperature  increase.  The  equation  for  a  time  constant,  t,  is 


T 


P  c 


L 


and  gives  the  time  to  reach  equilibrium 


(2-25) 


.  1 
T  -  T .  e 

r  1 


0.368 


(2-26) 


in  the  absence  of  radiation.  In  Fig.  2-8  the  generalized  time 
constant  for  an  arbitrary  set  of  conditions,  and  a  typical  material, 
MgF2,  0.07  in.  thick,  is  used  to  get  a  specific  time  constant.  For 
1  mm  thick  diamond  domes,  the  heat  transfer  coefficient  at  the 
stagnation  point  leads  to  very  small  Biot  nxmbers,  which  in  turn 
suggests  thermal  time  constants  of  the  order  of  1  sec  at  very  high 
altitudes.  In  these  conditions,  for  flight  times  of  about  10 
seconds  or  more,  diamond  missile  domes  will  indeed  reach  the 
recovery  temperature [ 4 ] . 

C.  AEROSTRUCTDRE  OF  THE  INFRARED  DOME 

As  will  be  discussed  in  later  sections,  the  self  emission  of 
the  dome  window,  due  to  the  high  temperatures  induced  by  aerokinetic 
heating,  is  an  aero-optical  effect  that  will  have  a  major  influence 
on  the  sensor  performance  in  a  high  speed  flight  situation.  Dome 
emission  is  going  to  depend  on  the  emissivity  of  the  dome  material 
among  other  variables.  An  accurate  prediction  of  the  dome  emissivity 
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should  consider  the  thickness  of  the  dome.  Therefore  a  discussion  of 
the  calculation  of  IR  dome  thickness  is  relevant.  The  dome's  ability 
to  withstand  the  thermal  shock  induced  by  transient  heating  on  a 
fly-out  trajectory  is  one  of  the  major  problems  encountered  in  the 
design  of  high  supersonic  IR  seeker  domes.  The  subject  of  thermal 
shock  resistance  of  IR  domes  has  been  very  well  characterized  and 
analysed  by  Dr.  Claude  Klein.  The  following  sub-sections  summarize 
some  of  his  published  work  on  the  subject [4] [10] [25] . 

1 .  Thermal  Shock  Resistance  of  the  IR  Dome 

One  of  the  most  constraining  factors  in  the  design  of  IR 
missile  window  domes  is  the  one  attributed  to  brittle  fracture 
caused  by  thermal  shock [10] ;  the  fracture  is  induced  by  tensile 
stresses  on  the  "cold"  side  corresponding  to  the  inner  surface  of 
the  dome.  These  stresses  are  caused  by  the  sudden  temperature 
gradients  generated  by  aerodynamic  heating,  which  in  conjunction 
with  thermal  expansion,  produce  compressions  on  the  outer  surface 
and  tensions  on  the  inside. 

Starting  with  the  expression  for  the  maximTom  stress  experienced 
by  a  clamped  plate,  or  a  complete  sphere. 


a 


max 


2 


a.E 

(1-v) 


AT  , 


(2-27) 


if  there  is  a  linear  temperature  variation  across  the  thickness  and 
if  both  are  free  to  expand: 


AT  -  |r„„-r  I 

'  OW  w' 


(2-28) 


where  AT  refers  to  the  temperature  difference  between  the  outer  and 
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where  AT  refers  to  the  temperature  difference  between  the  outer  and 
inner  surfaces.  In  a  transient  situation,  AT  is  a  strong  function  of 
time  and  depends  critically  upon  the  heat-transport  process  as 
characterized  by  the  Biot  number,  Bi .  If  the  condition  Bi<l  applies, 
AT  peaks  when  the  backface  temperature  starts  to  rise,  which  occurs 
long  before  the  front  reaches  the  stagnation  temperature,  and,  thus, 
leads  to 


AI  -  (T  ,-T.  )  (Bi) 

P  St  IW'  '  '  St 


(2-29) 


In  a  first  approximation,  therefore,  a  thermally -thin  dome  will 
be  subjected  to  a  peak  stress  that  amounts  to  [10] 


a£ 


max'p 


2  (1  -  V) 


(Si)  ,  (T  ,-r.  ) 

'  St  '  St 


(2-30) 


where  refers  to  the  Biot  number  at  the  stagnation  point,  T^^ 
is  the  relevant  stagnation  temperature,  and  must  be  set  equal  to 
the  wall  temperature  at  the  onset  of  the  shock.  Assessing  the 
thermal  shock  resistance  (TSR)  then  boils  down  to  comparing  |a^|p 
with  the  allowable  flexural  stress,  a^/i,  which  is  best  done  by 
defining  TSR  as  follows : 


TSR 


10  log 


a^/4 


a 


max 


(2-31) 


if  the  results  are  to  be  quantified.  The  stress  ratio  can  be 
expressed  as  a  product  of  three  terms. 
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a^/4  k 


The  first  term  involves  intrinsic  material  properties  only,  the 
second  term  depends  upon  both  the  diameter  of  the  dome (through  the 
heat  transfer  coefficient)  and  its  thickness,  while  the  third  term 
reflects  the  severity  of  the  aerothermal  environment. 


Figure  2-9.  Schematic  of  a  dome 
configuration[ 10], 

2 .  Pressure  Induced  Stresses 

At  Mach  numbers  M^3 ,  the  pressure  distribution  on  a  spheroidal 
shell  varies  essentially  as  the  cosine  square  of  the  stagnation 
9^1151®,  which  implies  that,  for  worst— case  type  calculations,  we  may 
assume  that  the  aerodynamic  pressure  on  a  truncated  hemisphere  is 
uniform  and  equal  to  the  stagnation  pressure  p^^..  In  high 
supersonic  flight,  a  fully  vented  IR  dome  as  in  Fig.  2-9,  must, 
therefore,  be  able  to  withstand  a  pressure  load [10] 
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Ap  -  P^j-P.  -  P^[  (p^^/pJ-1]  , 


(2-33) 


where  p^JP^  refers  to  the  pressure  ratio  across  the  bow  shock  and 
along  the  stagnation  streamline: 
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The  external  pressures  give  rise  to  a  complex  stress  pattern 
that  reflects  the  nature  of  the  edge  attachment,  i.e.,  the 
stiffness  of  the  dome  interface  with  the  support  ring.  Two  limiting 
cases  will  be  considered [ 10 ] ;  (a)  the  simply  supported  shell  that  is 
free  to  bend  when  subjected  to  the  action  of  a  uniform  normal 
pressure,  and  (b)  the  shell  with  built-in  (or  clamped)  edges  that  is 
restrained  from  moving  in  the  "radial"  direction.  For  simply 
supported  shells,  the  maximum  stress  induced  by  bending,  is 

tensile  and  obeys  the  equation 
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—  sin (Q) 

cos  (Q)  -  1 

[1  N 

^  J 

(2-35) 


provided  the  dome  aspect  ratio,  L/R,  satisfies  the  condition 


sin^(6)  L  sin^ (Q) 

-  2  — £  -  , 

12  R  1.2 


(2-36) 


where  6  designates  the  hemisphere  truncation  angle  (see  Fig. 2-9) . 
For  edge-constrained  shells,  the  predominant  stress  is  compressive 
and  localized  at  the  base  of  the  dome,  near  the  interface.  It 
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relates  to  the  aspect  ratio  of  the  dome  in  a  very  simple  relation 


sin^ (Q)  L  sin^(6) 

I2  R  3 


if  the  condition 


(2-37) 
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is  met.  Evaluating  Eg. (2-35)  for  several  angles  6  shows  that,  on  a 
relative  scale,  the  bending  stress  depends  solely  on  the  wall 
thickness  (the  stress  increases  as  the  thickness  decreases)  and  can 
be  approximated  by  means  of  the  relation 


a 

max 

Ap 


0.581 


(2-39) 


as  long  as  the  conditions  (2-37)  are  satisfied.  It  is  also  seen  that 
this  stress  represents  a  worst  case  in  the  sense  that  constraining 
the  edge  displacement  would  decrease  the  stresses,  especially  for 
low  aspect  ratio  domes.  Because  of  the  thermal  expansion  mismatch  at 
the  interface  of  the  infrared  dome  and  its  support  structure,  a 
"floating"  attachment  scheme  that  decouples  the  dome  from  the 
missile  body  would  be  desirable,  and  such  optimized  attachment 
should  result  in  pressure-induced  tensile  stresses  as  given  by 
Eg. (2-39) .  The  maximum  tensile  stresses  failure  criterion,  i.e., 

^  (2-40) 
where  refers  to  the  nominal  flexure  strength  of  the  dome  material. 
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and  SF  is  an  appropriate  safety  factor,  then  yields  the  minimum 
thickness  required  to  prevent  fracture.  Since  currently  available  IR 
transmitting  dome  materials  are  brittle  ceramics,  they  may  exhibit 
a  wide  range  of  fracture  probabilities  under  apparently  identical 
loadings,  which  suggests  safety  factors  of  at  least  four  (4)  if 
is  the  characteristic  flexural  strength  derived  from  Weibull 
statistical  analysis.  It  is  concluded  that  the  dome  thickness 
required  for  withstanding  a  pressure  load  Ap  with  good  probability 
of  survival  can  be  obtained  from  the  following  equation: 


L  . 

min 

R 


1. 75 


2/3 


(2-41) 


D.  RADIATION  PROPERTIES  OF  SURFACES 

In  analysing  the  interchange  of  radiant  energy  among  surfaces, 
it  is  necessary  to  know  their  emission,  reflection,  and  absorption 
characteristics .  Besides  the  obvious  need  for  transparency,  a 
critical  reguirement  for  an  IR  window  in  a  hot  environment  is  that 
it  emits  little  radiation  that  would  obscure  the  scene  being  viewed. 
In  the  next  sections  the  emission  properties  of  windows  will  be 
explored. 

1 o  Radi ome try 

In  describing  the  directional  transport  of  energy,  it  is 
^o^venient  to  deal  with  the  radiant  energy  (in  joules)  contained 
within  a  small  solid  angle  d<o.  This  naturally  leads  to  the  concept 
of  the  radiant  intensity.  The  radiant  flux,  0  (in  watts)  emitted 
per  unit  solid  angle  (cc)  by  a  point  source  in  a  given  direction  is 
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called  radiant  intensity,  I,  and  is  given  by [19] 


I  -  -  [  IV/sr]  . 

dco 


(2-42) 


Another  useful  radiometric  term  is  the  radiant  flux  density  at  a 
surface,  measured  in  watts  per  square  meter.  It  may  be  either 
emitted  from  a  surface  ,  in  which  case  it  is  called  radiant 
exitance,  M,  or  incident  onto  a  surface,  in  which  case  is  called 
irradiance,  E. 

The  inverse-square  law  of  radiation  from  a  point  source  is 
derived  by  calculating  the  irradiance  of  a  point  source  on  a 
spherical  surface  surrounding  the  point,  of  solid  angle  An  sr  and 
surface  area  4nr^ .  Hence 


il>  4ni  I 
A  4nr^ 


(2-43) 


The  radiance,  L,  represents  the  radiant  energy  per  unit  time 
and  unit  surface  area  that  passes  through  a  solid  angle  do  inclined 
at  an  angle  6  with  respect  to  the  surface  normal.  This  situation  is 
represented  in  Fig.  2-10.  The  radiance  is  given  by 


(dA  cosQ) 


(2-44) 


In  the  case  of  spectral  quantities (wavelength-dependent 
distribution) ,  a  subscript  A  is  affixed  to  the  radiometric 
quantities,  whereas  for  total  radiation  these  subscripts  are 
omitted.  In  general,  the  direction  of  the  intensity  is  specified  by 
angles  0  and  0  of  a  spherical  surface  as  shown  in  Fig.  2-10.  In  the 


case  of  a  Lambertian  surface,  the  intensity  for  a  certain  radius  of 
the  surface  varies  in  accordance  with  the  following  relation. 


I C©  j  -  1(0)  cos  6  . 

In  this  case  the  radiance  is  independent  of  the  viewing  angle,  and 
the  surface  is  said  to  be  perfectly  diffuse.  The  intensity 
distributions  are  related  to  the  corresponding  hemispherical 
radiation  fluxes  by  integration  over  the  hemisphere  surrounding  the 
surface.  As  an  example,  for  a  flat  surface,  the  radiant  power  0 
leaving  the  surface  is  connected  with  the  radiant  intensity  at  the 
surface  as  follows: 


(2-45) 


<i>  .  1(0)  cos  0  dco  [W] 

hem 


(2-46) 


in  which  the  subscript  'hem'  denotes  the  hemisphere.  Similar 
relationships  apply  for  the  incident  and  reflected  radiant  fluxes. 


ft 


Figure  2-10.  Angle  Notation. 
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2 .  Blackbody  Radiation 

The  ideal  blackbody  is  a  perfectly  efficient  and  diffuse 
emitter  of  radiant  energy.  No  body  at  the  same  temperature  can  emit 
more  radiation  at  any  wavelength  or  into  any  direction  than  can  a 
blackbody  at  the  same  temperature.  It  follows  that  a  blackbody  is 
also  a  perfect  absorber:  all  radiation  incident  on  a  blackbody  will 
be  absorbed  regardless  of  wavelength  or  angle  of  incidence.  The 
spectral  radiant  exitance  of  a  blackbody  is  given  by  the  Planck 
equation, 


2nhc^ 

A" 


I 

g  hc/XkT  2 


(2-47) 


where  the  physical  constants  h,  c,  and  k  represent  the  Planck 
constant,  the  speed  of  light,  and  the  Boltzmann  constant 
respectively.  After  substituting  the  values  of  these  constants  in 
the  Planck  equation,  the  result  is 


^x- 


3.  745  X  10^ 
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(2-48) 


where  A  is  in  micrometers  and  T  is  in  kelvins .  The  spectral  radiant 
exitance  for  three  background  temperatures  is  plotted  in  Fig.  2-11. 

The  peak  exitance,  shifts  toward  shorter  wavelengths  with 

increasing  temperature  according  to  the  Wien  displacement  law, 
expressed  as 

X  T  .  —  .  2.88  X  10^  [p/n-K]  .  '  (2-49) 

max  5j[,  L  K  J 

Integration  of  the  spectral  radiant  exitance  of  Eq. (2-47)  over  all 
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wavelengths  gives  the  total  radiant  exitance  given  by 


M  .  OT^ 


(2-50) 


known  as  the  Stefan-Boltzmann  law,  with  a  as  the  Stefan-Boltzmann 
constant,  equal  to  5 . 67xl0‘®W/m^-K^ . 

The  emissivity  e  expresses  quantitatively  the  fact  that  real 
surfaces  always  emit  less  radiation  then  a  blackbody.  The  blackbody 
emissivity  is  expressed  as  the  ratio  of  the  total  radiant  exitance 
of  the  body  to  the  total  radiant  exitance  of  a  blackbody  at  the  same 
temperature, 


e  (Tj  . 


(2-51) 
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Figure  2-11.  Blackbody  radiation  spectral  distribution  at  three  background 
temperatures  [20 ]. 

3 .  Directional  Emissivity 

The  angular  emissivity  is  defined  as  the  ratio  of  the  radiance 
in  a  particular  direction  to  the  radiance  of  a  blackbody  in  that 
same  direction.  In  general,  for  arbitrary  surface  conditions,  the 
distribution  of  the  emitted  radiance  L  depends  on  two  angles  6  and 
(j)  as  shown  in  Fig.  2-10.  On  the  other  hand,  the  blackbody  radiance 
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Lbb  is  uniform.  Consequently,  one  can  write 


£  (T,e,(p)  .  Jl.  .  s  , 

Ij  ryry 


(2-52) 


Eq.{2-52)  states  that,  for  a  diffuse  surface,  the  directional 
emissivity  is  independent  of  angle  and  is  equal  to  the  hemispherical 
emissivity . 

The  directional  emissivity  characteristics  of  real  surfaces  are 
revealed  by  measurements  performed  on  a  variety  of  materials  by  E. 
Schmidt  and  E.  Eckert [11].  One  of  their  measured  distributions  is 
reproduced  in  Fig.  2-11.  The  figures  are  in  the  form  of  polar 
diagrams  in  which  the  concentric  circles  are  contours  of  constant 
emissivity  e  and  the  radial  lines  are  contours  of  constant  angle  6. 
The  surfaces  employed  in  the  experiments  were  presumably  isotropic, 
and  no  dependence  on  the  angle  $  is  indicated.  The  temperatures  of 
the  various  test  surfaces  were  typically  a  few  hundred  degrees 
Fahrenheit.  The  directional  emissivity  results  were  determined  from 
total  radiation  measurements . 

Figure  2-11  shows  that,  for  nonmetals,  the  directional 
emissivity  remains  essentially  uniform  for  inclination  angles 
between  0°  and  50-60°;  then  it  drops  off  quite  rapidly  to  zero.  In 
the  region  of  near  uniformity,  the  emissivity  values  are  in  the 
range  from  0.8  to  0.95,  which  is  typical  for  dielectrics. 
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Figure  2-12.  Distribution  of  the  total  directional  emittance  for  several  dielectrics[l  1 ]. 

The  directional  emissivity  characteristics  displayed  by 
metallic  surfaces  are  somewhat  different  from  those  just  discussed. 
In  the  range  of  angles  between  0°  and  3  0-40°,  s  remains  quite 
uniform.  However,  at  larger  angular  inclinations  s  increases  quite 
sharply.  The  percentage  increase  is  more  pronounced  when  the  level 
of  emissivity  is  low.  At  angles  approaching  90°  the  directional 
emissivity  should  approach  zero . 

The  directional  absorptivity  is  the  fraction  of  the  radiant 
energy  incident  on  the  surface  from  a  specific  direction  that  is 
absorbed.  In  accordance  with  Kirchoff's  law,  the  directional 
monochromatic  absorptivity  and  the  directional  monochromatic 
emissivity  are  equal  provided  that  the  incident  beam  is  uniformly 
polarized.  That  is, 

(Q,<p)  .  (2-53) 
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4.  Optical  Properties  of  Dielectric  Window  Materials 
In  the  infrared,  there  are  two  main  contributors  to  the  complex 
index  of  refraction  of  insulating  ceramics.  They  are  electronic 
transitions  and  lattice  vibrations  or  phonons  {a  quantized  packet  of 
vibrational  energy)  [12]  .  The  electronic  transitions  contribute  to 
the  real  part  of  the  index  of  refraction  only  in  the  foinn  of  a 
frequency- independent  parameter  with  weak  temperature  dependence. 
Most  of  the  frequency  and  temperature  dependence  of  the  complex 
index  of  refraction  at  mid- infrared  frequencies  comes  from  phonons. 

Single-phonon  absorption  processes  are  due  to  the  excitation  of 
a  phonon  upon  the  absorption  of  a  photon.  These  processes  are  very 
strong,  since  the  photons  can  couple  directly  to  the  fundamental 
lattice  vibrations  of  a  material.  Multiphonon  absorption  processes 
are  due  the  excitation  of  more  than  one  phonon  upon  absorption  of  a 
photon.  These  processes  become  important  in  the  transparent  regions 
of  a  material  where  one-phonon  absorption  processes  do  not  dominate. 
Difference  bands  of  multiphonon  absorption  occur  in  the  millimeter 
spectral  region,  and  sum  bands  occur  in  the  infrared  region.  Sxom 
band  processes  typically  dominate  infrared  emissivity [12] . 

The  behavior  of  the  absorption  coefficient  and  refractive  index 
of  a  typical  optical  material  is  indicated  schematically  in  Fig.  2- 
12.  In  general,  materials  absorb  light  strongly  in  the  infrared  and 
ultraviolet  regions,  near  wavenumbers  of  10^  and  10^  cm"^, 
respectively,  in  the  upper  part  of  Fig.  2-13  [12]  .  Infrared 
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Figure  2-13.  Absorption  coefficient  and  refractive  index  of  a  typical  optical  window  material.  Adapted 
from  Thomas [9]. 

radiation  stimulates  vibrations  of  the  atoms  in  a  solid  at  the 
characteristic  Reststrahlen  frequencies.  Ultraviolet  radiation 
promotes  electrons  from  filled  bands  to  higher  energy  unfilled 
bands .  Between  the  two  strong  absorption  regions  is  a  window  in 
which  the  material  is  transparent.  The  infrared  edge  of  the  window 
is  called  the  laultiphonon  region  because  the  transitions  responsible 
for  the  weak  absorption  arise  from  multiple  excitations  of  one  or 
more  vibrational  modes (phonons) .  The  ultraviolet  edge  of  the  window 
is  called  the  Urbach  tail. 

The  optical  properties  of  non-magnetic  dielectric  window 
materials  can  be  specified  completely  by  the  complex  index  of 
refraction,  n(v,T)  or  complex  permittivity,  6'(v,T)  where  v  is  the 
wavenumber  (reciprocal  wavelength)  and  T  is  temperature [12 ] .  The 
relationship  between  n  and  e  is 
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n(V,T)  .  v'e  (V ,T)  .  n(V,T)  -  ik(V,T) 


(2-54) 


where 


&(V)  ■  (n^  -  k^)  -  2ink  -  -  ie^  ,  (2-55) 

n  is  the  index  of  refraction,  and  k  is  the  index  of  extinction. 

The  magnitude  of  the  reflectance  of  the  medixam  (at  normal 
incidence)  is 


2  1,2 


\R\  ■ 


k 
k^ 


(2-56) 


The  absorption  coefficient  is  defined  by 


^Jv^T)  .  4uk(v,T) 


The  radiance  from  a  material,  ('v,T),  is  expressed  as 


L^(V,T)  .  £  (Q,^,T)L^^(V,T) 


(2-57) 


(2-58) 


where  £(e,v,T)  is  the  spectral  emissivity  of  the  material,  Ljgs(v,T) 
is  the  radiance  of  a  blackbody,  T  is  the  temperature,  v  is  the 
spectral . frequency  and  9  is  the  angle  of  incidence  (i.e.  the  angle 
between  the  surface  normal  and  the  incident  ray) .  In  general,  the 
emissivity  is  composed  of  two  components:  a  specular  and  a  diffuse 
component.  For  our  purposes,  the  diffuse  component  is  negligible 
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since  the  materials  have  smooth  flat  surfaces  and  negligible  bulk 
scattering . 

Based  on  the  Stefan-Boltzmann  law  and  Kirchof  f '  s  law  of 
^^diation^  the  sgecular  component  of  emisszvit^  can  be  expressed 
as  [12] 


/7  D/Q  rn\  \  /-}  t/cosQ  . 

Eje.v.T)  .  _Ll-  (i  .  e  ‘ _ ^ 

1  - 


(2-59) 


where  R  is  the  single-surface  Fresnel  power  reflectance  for 
unpolarized  light,  and  t  is  the  window  thickness.  In  general  and 
R  have  angular  dependencies  consisting  of  sine  and  cosine  functions. 
Near  normal  incidence  (i.e.  near  9=0°  ),  the  angular  dependencies  of 
and  R  can  be  neglected,  since  the  sine  and  cosine  are 

approximately  constant. 

Usually,  a  measurement  of  \r\  is  not  sufficient  to  characterise 
e.  Therefore,  it  is  convenient  to  describe  the  form  of  the  complex 
permittivity.  The  classical  pole-fit  model  represents  the 
fundamental  lattice  vibrations  or  single-phonon  transitions  in  a 
material.  This  feature  is  important  in  representing  the  index  of 

refraction  at  all  optical  frequencies,  but  the  absorption 

coefficient  obtained  from  the  model  is  valid  only  in  the  vicinity  of 
the  one-phonon  frequencies.  The  absorption  coefficient  needs  to  be 
known  in  the  region  of  transparency  as  well,  requiring  a  description 
of  multiphonon  absorption  bands  in  which  more  than  one-phonon 

excitation  occurs  upon  the  absorption  of  a  single  photon. 
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E.  DOME  EMISSION  AND  EMISSIVITY  AT  HIGH  TEMPERATURES 


As  the  IR  dome  heats  up,  it  will  emit  graybody  radiation  over 
the  entire  electromagnetic  spectrum,  which  may  adversely  impact  the 
operation  of  the  seeker  system  in  a  high-speed  flight  situation 
[2]  [4]  [13]  .  It  is  recognized  that  saturation  of  the  detectors  by 
radiation  from  aerodynamically  heated  domes  is  not  a  major  issue  but 
the  added  background- induced  noise  will  lower  the  signal-to-noise 
ratio  for  detection,  decreasing  the  effective  performance  of 
the  system  in  terms  of  detection  range  as  well  as  tracking  ability. 
An  accurate  prediction  of  the  dome  emission  requires  a  good 
knowledge  of  the  dome  emissivity  properties  at  the  high  operating 
temperature.  The  emissivity  can  be  obtained  by  models  or  by 
experimental  data . 

The  effect  of  temperature  on  the  emissivity  properties  of 
several  dome  materials  is  shown  in  Fig.  2-14.  The  absorption  edge 
of  an  infrared  window  material  shifts  to  a  shorter  wavelength  as  the 
temperature  increases  (Fig.  2-14).  This  occurs  as  the  weak 
absorptions  present  in  this  region  begin  to  dominate. 

When  temperature  of  a  material  increases (see  Fig.  2-15),  so  does 
its  emissivity.  Fig.  2-16  shows  the  effect  predicted  for 
sapphire [9]  .  As  shown  in  Fig.  2-15  the  sapphire  window  at  1000°C 
(1273  K)  emits  7.4  times  as  much  radiant  energy  as  the  same  window 
at  500°C  (773  K)  if  the  emissivity  is  independent  of  T[9]  .  Near  the 
absorption  edge  the  emissivity  increases  with  temperature,  so  the 
actual  emission  is  even  more  than  7.4  times  the  lower  temperature 
value . 
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Wavelength  (^im) 


{d)  Germanium  (single  crystal,  1.14  mm  thick) 


Wavelength  (^m) 


{b)  Silicon  (single  crystal,  4.16  mm  thick) 


Wavelength  (/im) 


Figure  2-14.  Emissivity  of  several  dome 
materials  [2]. 

At  elevated  temperatures,  it  is  emission,  not  absorption,  that 
limits  the  performance  of  an  optical  window. 

Figure  2-17  shows  the  measured  high- temperature  spectral 
emissivity  of  a  sapphire  sample (Michael  Thomas,  APL[12]). 
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Figure  2-15.  Sapphire  absorption  coefficient  as  a  function  of  temperature  [9]. 


Figure  2-16.  Temperature  dependence  of  infrared  transmission  of  zinc  sulfide  and  sapphire[9] . 


Emissivity 


Figure  2-1 7.  Emissivity  of  0.20  cm  thick  sapphire  calculated  from 
measured  absorption[9] . 


Figure  2-18.  High-temperature  spectral  emissivity  of  a  hemex  grade  sapphire  o-ray  sample  (12. 7  mm 
diameter,  1. 0  mm  thick)  over  the  temperature  range  of 950  K  to  1930K.  [12]. 
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III.  GOVERNING  EQUATIONS  OF  FLUID  MECHANICS 

In  order  to  compute  compressible  viscous  flow  around  a  body, 
the  continuity,  momentum  and  energy  equations  must  be  solved 
simultaneously.  Here  the  vector  and  the  conservative-law  form  of 
the  compressible  Reynolds -averaged  Navier-Stokes  equation  is 
presented. 

A.  CONTINUITY  EQUATION 

The  continuity  equation  expresses  the  Conservation  of  Mass  law. 
This  law  applied  to  a  fluid  passing  through  an  infinitesimal  control 
volume  fixed  in  space  yields  the  following  equation  of  continuity 

dp 

—  .  V  ■  (  pv  )  .  0  .  (3.1) 


Here  p  is  the  fluid  density  and  H  is  the  fluid  velocity  vector. 
Eq.  (3-1)  states  that  the  net  mass  flux  through  a  control  volume 
surface  must  be  equal  to  the  time  rate  of  change  inside  the  control 
volume.  For  a  three-dimensional  Cartesian  system,  where  u,  v,  w 
represent  the  x  ,y,  z  components  of  the  velocity  vector,  Eq. (3-1) 
becomes 


dt 


d 

—  (pw)  -  0 
dz 


(3-2) 
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B .  MOMENTUM  EQUATION 

The  momentiom  equation  expresses  Newton's  second  law  as  applied 
to  a  fluid  element  passing  through  an  infinitesimal,  fixed  control 
volume  in  space.  The  momentum  equation  is: 


pv  )  ,  V  ■  pw  .  pf  .  V  ■  .  (3.3) 

The  first  term  in  Eq. (3-2)  represents  the  time  rate  of  change 
of  momentiim  per  unit  volume  in  the  control  volume.  The  second  term 
is  a  vector  representing  the  momentum  convection,  per  unit  volume, 
through  the  control  surface.  The  first  term  on  the  right-hand  side 
is  the  body  force  per  unit  volime  and  Hij  is  the  stress  tensor  given 
in  compact  tensor  notation  by 


n. 


-p5,^  •  n 


du  . 
1 

dx  . 

J 


du  . 

_ J 

dx  . 


du 

dx, 


( i  f  j  ^  ]c  -  1.  f  2  f  3  ) 


(3-4) 


where  dy  is  the  Kronecker  delta  symbol  (6^^  =  1  if  i  =  j  and  5^^  =  0 
if  i  J  )  ;  Uj  ,  U2  ,  U3  represent  the  three  components  of  the 
velocity  vector  V;  ]i  is  the  coefficient  of  viscosity  (dynamic 
viscosity)  and  X  is  the  second  coefficient  of  viscosity. 

By  substituting  Eq. (3-4)  into  Eq. (3-3)  and  writing  the  Cartesian 
component  velocities  in  the  directions  x,  y,  z  as  u,  v  and  w,  for 
clarity,  the  Navier-Stokes  equation  is  obtained.  In  component  form 
the  Navier-Stokes  equations  are : 
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(3-5) 


d  \ 

^  dv 

d 

iJ 

2—  - 

du 

dy 

dx  ^  1 

dx 

dy) 

dy 

[^1 

, 

dx 

~dz  J 

(3-6) 


dp  a 

Qf  ,  ^  ^  ^ 

*  az  dx 


dv  dw 


d  2  I  dw  du 

- m  2 -  -  - 

dz  3  \  dz  dx 


C .  ENERGY  EQUATION 

The  energy  equation  is  obtained  by  applying  the  First  Law  of 
Thermodynamics  to  a  fluid  passing  through  an  infinitesimal,  fixed 
control  volume 


do  tt 

- ^  .  V  •  V  .  —  .  V  •  g-  .  pf  •  V  .  V  •  (U.  ■  V) 

at  at 


(3-8) 


where  is  the  total  energy  per  unit  volume.  The  first  term  in  left 
hand-side  of  Eq. (3-8)  represents  the  rate  increase  of  total  energy 
per  unit  volume  in  the  control  volume  while  the  second  term 
represents  the  rate  of  total  energy  flux  per  unit  volume  through  the 
control  surface.  The  first  term  on  the  right  hand-side  is  the  rate 
of  heat  produced  by  external  agents  where  Q  is  the  heat  added  to 
the  system.  V-g  is  the  rate  of  heat  conduction  through  the  surface. 
Fourier's  law  for  heat  transfer  by  conduction  states  that  the  heat 
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transfer  q  can  be  expressed  as 

q  =  -kVr  (3-9) 

where  k  is  the  coefficient  of  thermal  conductivity  and  T  is  the 
temperature,  f  represents  the  body  forces  per  unit  mass.  pf-V  is 
the  work  done  on  the  control  volume  by  body  surfaces,  and  V -(11^^.  -V) 
is  the  work  done  on  the  control  volume  by  surface  forces . 

The  conservation  of  energy  for  a  three-dimensional  Cartesian 
coordinate  system  becomes 


dt 


-  -  P  (fxU  +  fyV  +  f^w)  +  (e^u  +  pu  -  UT^  -  vv^ 

^  ^  pv  -  ur^  -  vr^  ^  +  q^) 

^  (e,  ^wpw  -  -  vvy,  -WT,,  +  q^)  =  0 


where 


(3-10) 
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-f  . 

dw'^ 

dy\ 

(3-11) 
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D.  EQUATIONS  IN  NONDIMENSIONAL  VECTOR  FORM 

In  order  to  facilitate  the  numerical  implementation  the 
governing  fluid  dynamic  equations  are  written  for  a  body  fitted 
coordinate  system.  The  equations  are  transformed  to  a  generalized 
body  fitted  coordinate  system  using  the  following  transformations 

l(x,y,z),  r\  .  X\(x,y,z)  ,  ^(x,y,z)  (3.12) 


PHYSICAL  DOMAIN 


COMPUTATIONAL  DOMAIN 


Figure  3-L  Transformations  and  body  coordinate  system[15]. 


The  strong  conservation  law  is  retained  for  a  curvilinear 
coordinate  system  r)  ,  ^  and  the  governing  equations  become  [15] 
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d^q  .  B^iE  .  e\)  .  d^{F  ,  F^)  .  d^{G\  g\)  .  0 


(3-13) 


where 
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■ 
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(3-14) 


V  -  ri^.ri^u.ii^v.ri^w 
W  -  V*  I  w 


(3-15) 


where  U,  V,  and  W  are  contravariant  velocities  written  without 
metric  normalization. 

The  viscous  terms  are  given  by 
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(3-16) 
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Here  Cartesian  derivatives  are  expanded  in  terms  of  ^  ,  n  ,  ^ 
coordinates  by  the  chain  rule  relations  as 


(3-18) 


The  Cartesian  velocity  components  u,  v,  w  are  non- 
dimensionalized  with  respect  to  a„  (the  free-stream  speed  of  sound), 
density,  p,  is  referenced  to  p„  and  total  energy,  e,  to 
p^aj .  Pressure  is  given  by 


P  ■  (X-1)  [e^ 


X  p  (u^.v^.w^)  ] 


(3-19) 


and  throughout  y  is  the  ratio  of  specific  heats.  Also,  k  is  the 
coefficient  of  thermal  conductivity,  p  is  the  dynamic  viscosity, 
while  A  from  the  Stokes'  hypothesis  is  -(2/3 Ip. 
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The  metric  terms  and  the  Jacobian  of  the  transformation  are 
obtained  from  chain  rule  expansion  of  ,  etc.,  and  solved  for 

/  etc.  These  terms  are  given  by. 


■  J(\yz-yx^^^^ 


rix  ■  J(z^y^-y^^^) 

r]y-  J(x^z^-x^z^) 

T\^  .  J(y^x^.x^y^) 

■  -^T^x-^T^y-^^z 

■  -^x^x-yt^y-^T^z 


and 


r-l 


(3-20) 


(3-21) 


E.  THIN  LAYER  APPROXIMATION 

In  order  to  accurately  compute  the  normal  gradients  in  the 
boundary  layer,  it  is  necessary  to  make  the  normal  grid  spacing  very 
fine  close  to  the  body  surface.  In  the  streamwise  direction  the  flow 
gradients  are  not  large,  when  compared  to  those  in  the  normal 
direction,  and  the  resolution  of  the  grid  can  be  lower.  This  results 
in  a  mesh  where  the  cells  have  a  high  aspect  ratio  close  to  the  body 
and  a  smaller  aspect  ratio  when  approaching  the  outer  boundary  of 
the  computational  domain.  This  procedure  is  going  to  save 
considerable  computer  storage  and  computation  time  without 
sacrificing  the  accuracy  of  the  computation. 

In  a  thin-layer  approximation  all  viscous  derivatives  in  the  ^ 
and  ri  direction  (along  the  body)  are  neglected,  while  terms  in 
normal  direction  ^  are  retained,  (see  Fig. 3-1) .  Eq. (3-1)  thus 
simplifies  to 
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9^g  .  8/  .  d^F.  dG  .  Re-^8^  s' 


(3-22) 


where 


s'. 


J-^ 


0 


P  (C  C*  O  ^X 

P  r^". o  ^c  ^^x“r  ^y^r  ^y 
p  rd-  ^y-  ll)  rp/3;  rc,u^.  C3,v^.  iw^) 

(^i-  Cy-  ^l)  [0 . 5f2  (u^.v^.w^)  ^.KPr'^  (y-l)'^  (a^)  ^] 

.(p/3)  ((p/3)  (^^u^.(,^v^.^^w^) 


(3-23) 


Unlike  the  boundary- layer  model  the  thin  layer  accounts  for 
pressure  variations  in  the  normal  direction  and  does  not  require 
matching  an  inviscid  layer  with  a  viscous  layer. 

The  eddy  viscosities  were  computed  using  the  algebraic  Baldwin- 
Lomax  eddy  viscosity  model . 
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IV.  FLOW  FIELD  COMPUTATION 


The  flow  field  was  computated  in  order  to  obtain  the  temperature 

variation  on  the  surface  of  the  dome,  corresponding  to  the  IR 
window.  The  aerodynamic  heating  of  the  infrared  window  is 
responsible  for  a  degradation  of  the  seeker  infrared  sensor 
performance  at  supersonic  speeds . 


A.  METHODOLOGY 

For  the  space  discretization  the  Steger-Warming  flux  vector 
splitting  is  used  along  the  streamwise  direction.  Along  the 
circumferential  and  normal  directions  central  differences  with  added 
numerical  dissipation  are  used.  The  integration  is  performed  with  a 
two-factor  partially  flux  split  implicit  algorithm.  A  combination  of 
second  and  fourth  order  differences  have  been  used  for  smoothing 
along  the  central-differencing  directions.  Fourth-order  dissipation 
terms  are  added  explicitly  to  the  equations  to  damp  high  frequency 
growth  and  thus  serve  to  control  nonlinear  instability.  The 
resulting  implicit  code  is  second-order  accurate  in  space. 

A  body-fitted  grid,  clustered  in  the  radial  direction  was 
employed  to  resolve  the  viscous  layer  adjacent  to  the  body  surface. 
This  grid  was  obtained  by  means  of  a  hyperbolic  grid  generator. 

The  flow  field  was  computed  using  the  OVERFLOW  Version  1.6t 
code  developed  at  NASA  Ames  Research  Center  running  in  a  STARDENT 
workstation  installed  at  the  Computational  Fluid  Dynamics  Laboratory 
(Naval  Postgraduate  School) .  After  a  solution  was  obtained,  the 
original  grid  points  were  redistributed  (adapted)  based  on  the 
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initial  flowfield  solution,  and  a  nsw  solution  was  computad  using 
the  adapted  grid.  This  procedure  increases  the  accuracy  of  the  flow 
prediction  by  decreasing  the  effects  of  flowfield  discretization. 
This  procedure  is  particularly  important,  in  supersonic  flow 
computations,  where  the  existence  of  strong  gradients  and  shocks 
require  careful  distribution  of  grid  points  to  minimize  errors  and 
produce  more  accurate  flow  field  predictions . 

B.  BODY  GEOMETRY  AND  GRID  GENERATION 

The  flowfield  was  computed  around  three  distinct  body 
geometries,  namely  two  nose  shapes  (cases  I  &  II)  and  a  test 
geometry  (case  III)  .  The  case  III  geometry  had  the  shape  of  a 
hemisphere-cylinder.  The  computed  shock  stand-off  distances  and 
surface  pressure  distribution  of  this  latter  configuration  was 
compared  with  experimental  data  to  give  confidence  in  the 
computations . 

Several  C-0  type  grids  were  generated,  corresponding  to  the 
hemispherical-ogive  (case  I)  and  to  the  ellipsoidal-ogive  (case  II) 
forebody.  Longitudinal  cross  sections  of  the  computational  grids  for 
the  case  I  and  case  II  are  shown  in  Fig.  4-3  and  Fig.  4-4 
respectively.  The  case  I  body  consists  of  a  truncated  hemisphere 
forebody (1 . 9024  cm  long)  attached  to  a  3:1  tangent  ogive  (12.728  cm 
long)  attached  to  a  4  caliber  tangent  cylindrical  afterbody.  The 
case  II  body  consists  of  a  truncated  ellipsoidal  forebody (2 . 6943  cm 
long)  attached  to  a  3:1  tangent  ogive  (12.72  8  cm  long) 
attached  to  a  4  caliber  tangent  cylindrical  afterbody.  The 
hemispherical  and  ellipsoidal  portions  of  the  nose  correspond  to  the 


52 


IR  window  locations  respectively.  The  origin  of  the  physical 
coordinates  is  at  the  nose  tip  for  both  cases.  The  three- 
dimensional  body  fitted  grids  were  generated  using  the  HYPGEN  code. 
After  a  flow  solution  was  found  an  adaptative  grid  code [17],  SAGE, 
based  on  the  algebraic,  solution-adaptive  scheme  of  Nakahashi  and 
Deiwert,  was  used,  in  order  to  enhance  the  resolution  for  better 
capture  of  the  shock  front.  A  third  inviscid  grid  was  designed  to 
compute  the  flow  around  a  hemisphere-cylinder (case  III).- 
corresponding  to  the  test  case.  The  geometrical  shapes  of  the  bodies 
considered  are  as  follows, 

1.  Case  I.  Hemisphere-Ogive-Cylinder  (Fig.  4-1) 

This  nose  configuration  corresponding  to  the  hemispherical  dome 
is  defined  using  the  following  relation 

r  .  y'3.3166^-  (z-3.3166)^  ,  O^zsl.9024 

2  (z.3.3697)^ 

r.  — (z.3.3697)- - — - —  ,  1 . 9024szs  14 . 6303  (4-1) 

3  54 

r  .  6.0000  ,  ZS14.6303 

The  origin  of  the  coordinate  system  is  located  at  the  tip  of 
the  nose.  Fig.  4-3  shows  the  computational  grid  used  for  the 
calculations.  The  resolution  of  the  grid  is  increased  around  the 
shock  wave  where  the  flow  gradients  are  large. 
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2.  Case  II.  Ellipsoid-Ogive-Cylinder 

This  nose  configuration  corresponding  to  the  ellipsoidal  dome 
is  defined  using  the  following  relations: 
elipsoid  excentricity  =  0.9 


2-10.1376)^ 

r  . 

4.489, 

( - T 

f  0:^  2 . 6943 

(10.1376)^ 

2  , 

(2.2.577)2 

r  - 

3 

577  - 

54 

,  2. 6943£2il5 .4223 

r  - 

6.0000 

r 

22: 15 . 4223 

(4-2) 


Figure  4-4  shows  the  77x93x31  computational  grid  used  for  the 
calculations. 

3.  Case  III.  Hemisphere -Cylinder  (Fig.  4-2) 

This  body  configuration  corresponds  to  the  test  case.  Its  shape 
is  defined  by  the  following  relations: 


r  -  \/36.00^-  {z-6.00)^ 
r  .  6 . 00  , 


OszsG .  00 
ZkS.OO 


(4-3) 


C .  GRID  ADAPTATION 

Grid  adaptation  procedures  in  general  either  redistribute 
points  or  refine  distribution  (by  adding  points)  to  reduce  solution 
errors  by  minimizing  grid  discretization  errors.  Gnoffo[17]  first 
introduced  a  method  for  the  redistribution  of  grid  points  based  on 
local  flow  gradients .  This  method  is  analogous  to  finding  the 
equilibrium  position  of  a  system  of  springs  that  connect  each  node 
with  adjacent  nodes  and  whose  spring  forces  (or  tension  constant) 
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Figure  4-4.  Three  dimensional  Hemisphere-Ogive-Cylinder  computational  mesh. 
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is  proportional  to  the  local  error  or  weight  function.  The 
multidimensional  problem  of  grid  adaptation  is  split  into  a  series 
of  one-dimensional  (1-D)  problems  along  the  computation  coordinate 
lines.  The  reduced  1-D  problem  then  requires  a  tridiagonal  solver 
to  find  the  location  of  grid  points  along  a  given  coordinate  line. 
Multidirectional  adaptation  is  achieved  by  the  sequential 
application  of  the  method  in  each  coordinate  direction  (Fig.  4-5)  . 


initial  grid 


'  i 

:z 

:  [ 

— I- 

r~ . 

.  f 

f 

\ 

1 

marching  direction 
first  adaption 


Figure  4-5.  Splitting  of  2-D  adaptation  into  two  1 -directional  adaptations [17], 


In  the  adaptive  grid  procedure  outlined  by  Nakahashi  and  Deiwert 
(1985)  the  tension  force  directs  the  redistribution  of  points  to  the 
strong  gradient  regions .  To  maintain  smoothness  and  a  measure  of 
orthogonality  of  grid  lines,  torsional  forces  are  introduced  that 
relate  information  between  the  family  of  lines  adjacent  to  one 
another [16] .  The  smoothness  and  orthogonality  constraints  are 
direction-dependent,  since  they  relate  only  the  coordinate  lines 
that  are  being  adapted  to  the  neighboring  lines  that  have  already 
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been  adapted.  This  implies  that  the  solutions  are  nonunique  and 
depend  on  the  order  and  direction  of  adaptation.  Nonunigueness  of 
the  adapted  grid  is  acceptable  since  it  makes  possible  an  overall 
and  local  error  reduction  through  grid  redistribution. 

D.  RESULTS  AND  DISCUSSION 

Table  4-1  shows  the  computed  values  of  absolute  temperature, 
normalized  temperature  and  coefficient  of  pressure  along  the 
stagnation  streamline  obtained  by  the  solution  of  the  thin-layer 
Navier-Stokes  equations  about  the  hemisphere-ogive-cylinder  nose 
(Case  I)  at  zero  incidence  at  Mach  4.  A  77x7x35  grid  point  inviscid 
grid  was  used.  The  solution  was  considered  converged  after  1500 
steps.  Figs. 4-6  and  Fig. 4-7  show  the  plot  of  these  flow  quantities 
versus  distance  from  the  stagnation  point.  The  bow  shock  stand-off 
distance  taken  from  the  computed  flow  field  is  approximately  0.36cm, 
corresponding  to  a  A/R=.ll.  From  the  plot  of  Fig.  4-7  we  can  see 
that  the  temperature  jimp  across  the  shock  is  4.15  compared  with  the 
value  of  4.04  obtained  from  normal  shock  tables . 

Table  4-2  shows  the  computed  values  of  absolute  temperature, 
normalized  temperature  and  coefficient  of  pressure  on  the  surface  of 
the  nose  along  a  streamline  for  the  hemisphere-ogive  at  Mach  4  for 
the  same  conditions  of  Table  4-1.  Fig.  4-8  shows  the  computed 
temperature  profile  of  the  dome  surface  temperature.  The  computed 
normalized  stagnation  temperature  is  4.6  comparing  with  4.2  obtained 
using  Eq.  (2-7)  .  Fig.  4-9  shows  the  computed  dome  pressure 
distribution.  Fig.  4-10  shows  the  OVERFLOW  input  file  used  to 
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generate  the  values  of  Tables  4-1  and  4-2.  Fig.  4-11  shows  the 
three-dimensional  flow-field  around  the  missile  afterbody.  The  IR 
dome  corresponds  to  the  white  portion  of  the  nose  tip.  Fig.  4-12 
shows  the  density  field  in  front  of  the  IR  dome. 

The  flow  field  has  cylindrical  symmetry  at  zero  angle-of- 
attack.  With  missile  incidence  the  flow  is  no  longer  axisymmetric . 
Figs.  4-13  and  4-14  show  the  Mach  number  distribution  and  the  dome 
temperature  field  in  the  dome  for  an  angle  of  attack  of  5°.  Fig.  4-15 
gives  the  density  field.  It  was  obtained  using  the  viscous  adapted 
grid  of  Fig. 4-5. 

In  the  computation  of  the  flow  it  was  assumed  that  the  dome  had 
reached  thermal  equilibrium  and  the  wall  was  adiabatic.  In  order  to 
obtain  the  equilibrium  temperature  the  dome  should  be  divided  into 
elements  and  a  procedure  like  the  one  outlined  in  Eq. (2-24)  should 
be  carried  out  to  consider  radiation  and  conduction  losses .  The  heat 
transfer  coefficient  at  stagnation  is  given  by  Eq. (2-20) .  For  other 
elements  that  make  an  off-axis  angle  (angle  between  the  longitudinal 
axis  of  the  missile  and  the  line  between  the  dome  element  and  the 
dome  center  of  curvature) ,  6,  the  value  of  should  be  multiplied 
by  cos  (6)^''^.  In  this  study  the  dome  material  was  assumed  to  be 
sapphire  or  CVD (Chemical  Vapour  Deposition)  Diamond.  These  materials 
make  it  possible  to  fabricate  IR  domes  with  a  high  thermal  shock 
resistance (see  Appendix  A)  suitable  for  the  MWIR (Medium  Wave 
Infrared)  and  LWIR(Long  Wave  Infrared)  spectral  bands  respectively. 
In  Appendix  A  it  is  shown  that  for  a  diamond  window  at  Mach  4  a  dome 
thickness  of  1mm  is  enough  to  withstand  the  flight  conditions. 
Therefore  we  can  expect  that  conduction  losses  are  not  going  to  be 
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high  since  the  eirea  of  the  attachment  is  small.  Also,  since  the 
emissivity  of  these  materials  is  small  we  can  expect  small  radiation 
losses.  In  this  case,  considering  the  wall  adiabatic  might  be  a  good 
approximation. 


Table  4-1.  Computed  streamline  flow  field  properties  (Mach=4.0,  a=0).  \ 


z 

X 

y 

T(K) 

T/Tinf 

Cp 

O.OOOOOOOE+00 

0 . OOOOOOOE+00 

O.OOOOOOOE+OO 

1251.995 

4.333828 

1.769331 

-1.4072160E-02 

0 . OOOOOOOE+00 

O.OOOOOOOE+OO 

1251.995 

4.333828 

1.769331 

-2.845700lE’-02 

O.OOOOOOOE+00 

0 .OOOOOOOE+OO 

1244.208 

4.306874 

1,753680 

-4.3183811E-02 

O.OOOOOOOE+00 

O.OOOOOOOE+OO 

1237.092 

4.282242 

1 , 753231 

-5.3293991E-02 

O-OOOOOOOE+OO 

O.OOOOOOOE+OO 

1232.358 

4.265854 

1.735707 

-7.3827960E-02 

0 . OOOOOOOE+00 

O.OOOOOOOE+OO 

1231.002 

4.261161 

1.726179 

-3.9831688E-02 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

1228.993 

4.254206 

1.708787 

-0.1063549 

0 . OOOOOOOE+00 

O.OOOOOOOE+OO 

1227-959 

4.250627 

1-697652 

-0.1234523 

0 . OOOOOOOE+OO 

O.OOOOOOOE+OO 

1226.781 

4.246548 

1 . 588257 

-0.1411849 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

1226,759 

4.246473 

1-686133 

-0.1596204 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

1225.986 

4.243798 

1 . 679476 

-0.1788342 

0  - OOOOOOOE+OO 

O.OOOOOOOE+OO 

1225.059 

4.240590 

1.674541 

-0.1989112 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

1221.979 

4.229926 

1.657630 

-0.2199472 

0 . OOOOOOOE+OO 

0 .OOOOOOOE+OO 

1219.554 

4.221533 

1.653113 

-0.2420507 

0 . OOOOOOOE+OO 

O.OOOOOOOE+OO 

1212.904 

4 , 198514 

1.616354 

-0.2653450 

0 . OOOOOOOE+OO 

O.OOOOOOOE+OO 

1207.786 

4.180799 

1.510427 

-0.2899718 

0 . OOOOOOOE+OO 

O.OOOOOOOE+OO 

1177 . 408 

4.075641 

1.431967 

-0,3160940 

0 . OOOOOOOE+OO 

O.OOOOOOOE+OO 

997.7296 

3.453679 

0,7734069 

-0.3439003 

0  - OOOOOOOE+OO 

0 .  OOOOOOOE+OO 

696.3934 

2,410592 

0 . 3077108 

-0.3736109 

0 . OOOOOOOE+OO 

O.OOOOOOOE+OO 

473.6082 

1.539413 

0.1148101 

-0.4054836 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

337.6206 

1.168687 

2.9509827 

-0.4398234 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

291.7347 

1.009851 

2.3921142 

-0.4769938 

O.OOOOOOOE+OO 

0  .  OOOOOOOE+OO 

288.1387 

0.9974033 

-3.0908835 

-0.5174306 

0  .  OOOOOOOE+OO 

0  .OOOOOOOE+OO 

289.3298 

1.001526 

5.1899254 

-0.5616624 

0  .  OOOOOOOE+OO 

O.OOOOOOOE+OO 

288.7061 

0.9993671 

-2.2675097 

-0.6103359 

0 . OOOOOOOE+OO 

0  .OOOOOOOE+OO 

288.9403 

1 . 000178 

6.7792833 

-0. 6642512 

O.OOOOOOOE+OO 

0  .OOOOOOOE+OO 

288.3846 

0. 9999849 

-5 . 8636069 

-0.7244109 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

288.3838 

0.9999823 

-5.1110983 

-0-7920877 

0.  OOOOOOOE+OO 

0  .OOOOOOOE+OO 

288.3910 

1.000007 

2,2649765 

-0.3689228 

0  .  OOOOOOOE+OO 

O.OOOOOOOE+OO 

288.3875 

0.9999951 

-6.2584877 

-0.9570694 

O.OOOOOOOE+OO 

0  -  OOOOOOOE+OO 

288.3888 

0.9999998 

1.5391277 

-1.059409 

O.OOOOOOOE+OO 

0  .  OOOOOOOE+OO 

288.3887 

0.9999993 

-1 . 0430813 

-1.179889 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

238.3891 

1 .000001 

-1.4901161 

-1,324077 

O.OOOOOOOE+OO 

O.OOOOOOOE+OO 

288.3890 

1.000000 

-2.2351742 

-1.500000 

O-OOOOOOOE+OO 

0  .  OOOOOOOE+OO 

288.3888 

0. 9999998 

-1 . 4901161 
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Table  4-2.  Computed  hemispherical  dome  temperature  profile  (Mach=4.0,  a=0)  (adiabatic 
wall). 


2  X 

y 

T(K) 

T/Tratio 

Cp 

O.OOOOOOOE+OO  O.OOOOOOOE+00 

O.OOOOOOOE+OO 

1251.995 

4.333828 

1.769331 

9. 9657424E-05  -1 . 5062407E-02 

1.5062407E-02 

1267.223 

4.386541 

1.362640 

3.1045929E-04  -3 . 0966893E”02 

3.0966893E-02 

1299.271 

4.497478 

1.752715 

6.7767297E~04  -4 . 7760621E-02 

4.7760621E-02 

1280.968 

4.434122 

1.695181 

1.2542109E-03  -6 . 5493003E-02 

6.5493003E-02 

1249.678 

4.325809 

1.720564 

2.0872049E-03  -8 . 4215783E-02 

8.4215783E-02 

1257.582 

4.353167 

1.724007 

3.2189630E-03  -0.1039834 

0.1039834 

1249.013 

4,323505 

1.722371 

4.6833940E-03  -0.1248538 

0.1248538 

1246.343 

4 . 314264 

1.724451 

6.5123211E-03  -0.1468881 

0.1468881 

1244.572 

4 . 308134 

1.734758 

8.7461583E-03  -0.1701505 

0.1701505 

1242.742 

4.301799 

1.720822 

1.1447350E-02  -0.1947079 

0.1947079 

1241.205 

4.296480 

1.720025 

l,4705360E-02  -0.2206291 

0.2206291 

1239.578 

4.290847 

1.714795 

1.8604770E-02  -0.2479863 

0.2479863 

1237.847 

4.284853 

1.712322 

2.3207851E-02  -0.2768574 

0.2768574 

1236.551 

4.280367 

1.705234 

2.8605610E-02  -0.3073223 

0.3073223 

1233.810 

4.270882 

1.693714 

3.4932461E-02  -0.3394611 

0.3394611 

1230.839 

4-260595 

1.682505 

4.2309761E-02  -0.3733589 

0.3733589 

1225.294 

4.241401 

1.664191 

5.0861739E-02  -0.4091049 

0.4091049 

1219.603 

4.221704 

1.645974 

6.0756631E-02  -0.4467875 

0.4467875 

1213.211 

4.199575 

1 . 623081 

7.2170392E~02  -0.4864984 

0.4864984 

1207.618 

4.180215 

1.600384 

8.5286379E-02  -0.5283324 

0.5283324 

1202.201 

4.161465 

1.575645 

0.1003310  -0.5723819 

0.5723819 

1196,849 

4 .142937 

1.548760 

0.1175769  -0.6187352 

0.6187352 

1190.913 

4.122391 

1.517524 

0.1372850  -0.6674858 

0.6674858 

1184.456 

4 .100039 

1.482646 

0.1597211  -0.7187309 

0.7187309 

1177.089 

4.074539 

1.442449 

0.1852678  -0.7725453 

0.7725453 

1168.978 

4.046461 

1.397681 

0.2143953  -0.8289816 

0.8289816 

1160.507 

4 . 017141 

.1.350345 

0.2474633  -0.8881212 

0.8881212 

1151.573 

3.986216 

1.300280 

0.2845651  -0.9501302 
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Figure  4-7.  Plot  of  the  normalized  temperature  along  the  stagnation  streamline  for  the  computed  flow  jietd 
for  the  hemisphere-ogive-cylinder. 
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SGLOBAL 

CHIMRA-  .F. 
NQT  =  0 , 
$END 
$FLOINP 

ALPHA  =0.0 
SEND 
SVARGAM 
SEND 
SGRDNAM 

NAME  =  'Hem 

SEND 

SNITERS 

SEND 

SMETPRM 

IRHS  =  0, 
SEND 
STIMACU 

DT  =0.2 
SEND 
SSMOACU 

ISPECJ=  2, 
ISPECK=  2, 
ISPECL=  2, 
SMOO  =1.0 
EPSE  =0.3 
SEND 
SVISINP 

VISCJ  =  .F. 
NTURB  =  1 

ITTYP  =  1 

ITDIR  =  3 

JTLS  =  1 

JTLE  =  77 

KTLS  =  1 

KTLE  =  7 

LTLS  =  1 

LTLE  =  35 

TLPAR1=  0.3 
SEND 
SBCINP 

NBC  =  6 

IBTYP  =  15 

IBDIR  =  1 

JBCS  =  1 

JBCE  =  1 

KBCS  =  1 

KBCE  =  7 

LBCS  =  1 

LBCE  =  35 

SEND 
SSCEINP 
SEND 


NSTEPS=  600, 


FSMACH- 


RESTRT=  .T., 


=  6.00E6u 


NSAVE  =  50, 


520.000, 


Hemisphi-ogiv-  Cylinder  77x7x35  (A=0 , 0  deg,  M=4.00,  Re=6.00  mil.) 


ILHS  =  2, 


=  0.20,  ITIME=  1,  TFOSO  =  1.00, 


DIS2J 

DIS2K 

DIS2L 


DIS4  J 
DIS4K 
DIS4L 


0.02, 

0.02, 

0.02, 


VISCL  =  .T. , 


12, 

5, 

32, 

30 

-2, 

3, 

-3, 

-1 

1, 

1, 

1, 

77 

77, 

77, 

*77, 

77 

7, 

1  , 

1 

■7, 

7, 

7, 

7 

1  . 

1, 

35, 

1 

35, 

1 , 

35, 

35 

Figure  4-10,  Sample  OVERFLOW  Input  File. 
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Figure  4-12.  Gradient  of  index  "lens"created  by  the  compressed  gas  in  front  of  the  dome  at  Mach=4.0. 
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Figure  4~14.  Hemispherical  dome  temperature  field  for  a  Mach  4  seeker  with  5"  angle-of  attack  obtained 


by  solving  the  Thin-Layer  Navier-Stokes  equations. 
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E .  VALIDATION 

In  order  to  validate  the  computation,  a  comparison  between  the 
predicted  flow  field  about  the  missile  nose  and  experiments  was 
tried.  After  an  exhaustive  library  search  no  data  suitable  for 
comparison  with  the  computation  was  found  in  the  open  literature. 
Most  experiments  were  performed  at  hypersonic  speeds .  The  only 
suitable  experimental  data  to  compare  with  the  flow  field 
computation  were  ones  for  a  hemisphere-cylinder  at  M„  =  1.2  and  a 
Re  =  5.4  X  10®  ft.  For  this  comparison  a  61x7x41  inviscid  grid  was 
created.  The  experimental  data  was  taken  from  shadowgraphs  obtained 
in  the  AEDC  Aerodynamic  Wind  Tunnel  (IT) [18] . 

A  plot  of  the  shock  standoff  distance  at  both  the  body  axis  and 
the  wind  axis  for  a  hemisphere-cylinder  at  Mach  1.2  is  shown  in  Fig. 
4-17.  The  normalized  standoff  distance,  A/R  for  zero  incidence  is 
-1.26.  The  flow  was  computed  using  a  61x7x41  inviscid  grid.  The 
computed  normalized  standoff  distance  was  -1.0.  The  poor  result  is 
considered  to  be  due  to  lack  of  resolution  of  the  grid  near  the 
shock  front.  Due  to  restrictions  of  time  another  grid  was  not 
generated  to  try  another  comparison.  Fig.  4-16  shows  the  computed 
Mach  number  distributions  using  the  inviscid  Euler  solution. 
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Figure  4-16.  Hemisphere-Cylinder  Euler  solution  of  the  Mach  Number  flow  field  (61x7x41  grid). 


V.  INFRARED  SEEKER  DETECTION  PERFORMANCE  DEGRADATION 

The  major  effect  of  the  increase  in  dome  temperature  on  the  IR 
seeker  sensing  performance  at  high  supersonic  speeds  is  considered 
to  be  thermal  emission  to  the  detectors.  This  increased  dome 
radiation  contributes  to  the  noise  and  thus  increases  the  Noise 
Equivalent  Irradiance  (NEI)  to  the  detector  or  Noise  Equivalent 
Quantum  Flux  (NEQF)  to  the  detector  array,  decreasing  the  effective 
detectivity  of  the  vehicular  sensor  and  thus  decreasing  the 
detection  range  of  the  seeker. 

Initially  the  basic  figure  of  merit  parameters  will  be 
introduced.  Then  a  discussion  of  the  range  equations  for  non-imaging 
and  imaging  sensors  in  a  high  speed  flight  environment  will  follow. 
Several  computations  of  detection  ranges  for  a  staring  sensor  are 
carried  out  for  a  range  of  mean  dome  temperatures,  representing  a 
range  of  missile  speeds  up  to  Mach  4. 

A.  SENSOR  PERFORMANCE  PARAMETERS 

It  is  convenient  to  introduce  the  seeker  sensor  figures  of 
merit  before  starting  an  analysis  of  the  seeker  performance. 

1 ,  Basic  Figures  of  Merit 

A  basic  figure  of  merit  that  applies  to  all  detectors  with 
electrical  output  is  responsivity .  Responsivity  is  the  ratio  of  the 
output  (usually  in  amperes  or  volts)  to  the  radiant  input  (in 
watts) .  The  spectral  voltage  responsivity  of  a  detector  at  a  given 
wavelength  A  is  the  measured  voltage  output,  V^,  divided  by  the 
spectral  radiant  power  incident  on  the  detector,  <te(k)  [29]  : 
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(5-1) 


where  f  is  the  modulation  frequency. 

The  noise  equivalent  power  [NEP)  of  a  detector  is  the  required 
power  incident  on  the  detector  to  produce  an  average  signal  output 
equal  to  the  rms  noise  output.  Stated  another  way,  the  NEP  is  the 
signal  power  level  that  produces  a  signal-to-noise  ratio  of  1.  The 
voltage  signal  output  is 


V 

s 


31  $ 

V 


e 


so  the  signal-to-noise  ratio  is 


(5-2) 


S/N  . 


31  4> 

V  e 


(5-3) 


The  NEP  is  the  incident  radiant  power,  <l>g,  for  a  signal- tp-noise 
ratio  of  1 : 


3t  NEP 
1 .  — 
^RMS 

Solving  for  NEP  gives 


NEP  . 


V 

RMS 

V 


(5-4) 


(5-5) 


where  is  the  root -mean- square  noise  voltage  in  volts  and  is 

the  voltage  responsivity  in  volts  per  watt. 


76 


The  detectivity  D  of  a  detector  is  defined  as  the  reciprocal  of 
the  noise  equivalent  power: 


D  .  [P/  . 

NEP 


(5-6) 


A  more  useful  performance  parameter  is  the  normalized  detectivity  D* 
(pronounced  dee-star)  which  normalizes  D  to  the  detector  area  and 
noise  bandwidth,  Af,  as  follows, 


D-  .  D^A  Af 


NEP 


(5-7) 


given  in  units  of  cm-Hz^''^/W. 

The  NEP  can  be  divided  by  the  area  of  the  collecting  aperture,  A^, 
leading  to  a  parameter  called  noise  equivalent  irradiance  {NEI)  [27]  , 


NET  . 


NEP 


(5-8) 


2 .  Noise  Equivalent  Temperature  Difference  (NETD) 

The  most  widely  used  radiometric  measure  of  the  ability  of 
a  military  mid-IR  system  to  discriminate  small  signals  in  noise  is 
the  Noise-Equivalent  Temperature  Difference  or  Noise-Equivalent 
Differential  Temperature.  As  defined  by  J.  M.  Lloyd  [22]  ,  "  The  NETD 

is  the  blackbody  target-to-background  temperature  difference  in  a 
standard  test  pattern  which  produces  a  peak-signal  to  rms-noise 
ratio  (SNR)  of  one  at  the  output  of  a  reference  electronic  filter 
when  the  system  views  the  test  pattern". 

In  terms  of  system  parameters  the  ratio  of  the  differential 
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signal  voltage  to  the  detector  rms  noise  voltage  in  the  detector  can 
be  expressed,  for  a  rectangular  detector,  as [22] 


AV 

S 

V 


.  AT- 


a(3A^ 


n/. 


abAf„ 


/ 


dT 


D- (Xj  TJX)  dX 


(5-9) 


where  a,  0  =  detector  angular  subtenses  [radian] 

a,  b  =  detector  dimensions  [cm] 

D* (A)  =  specific  detectivity  as  a  function  of 
wavelength  ,  evaluated  at  the  electrical 
frequency  at  which  the  noise  voltage  spectrirm  of  the 
detector  is  normalized  to  unity  [cm  Hz^''^  /  W] 

Afi^  =  equivalent  noise  bandwidth  [Hz]  of  the  NETD  test 
reference  filter 

Aq  =  effective  area  of  the  collecting  aperture  [cm^] 

M^(A)  =  spectral  radiant  exitance  [W/pm-cm^] 

To  (A)  =  spectral  transmittance  of  the  optics. 

The  NETD  is  defined  as  the  temperature  difference  for  unity  SNR,  so 
we  can  set  AV^/V^  -  1  and  solve  for  the  resulting  AT,  which  is  by 
definition  the  NETD: 


NETD  ■ 
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n  ,Jab  _ 

d  M, 

- -  D-  (X)  T  (X)  dX 

dT  ° 


(5-10) 


The  spectral  radiant  exitance,  ,  is  defined  by  Planck's 

blackbody  radiation  law 
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where  Cj  =  3.7415  x  10^  [watt  ijm^  /  cm^  ] 

C2  =  1.4388  X  10^  [urn  K]  . 

The  noise  equivalent  temperature  difference,  NETD,  representing 
the  total  output  temporal  and  pixel (spatial)  noise  at  a  given 
background  temperature  converted  into  an  equivalent  temperature 
difference  at  the  background  scene,  is  a  very  effective  figure  of 
merit  for  a  staring  focal  plane  array  (FPA) . 

3  .  Noise  Eqpiivalent  Quantum  Flux  for  a  CCD  Array 
To  develop  such  an  array  figure  of  merit  (FOM) ,  an  appropriate 
performance  measure  for  a  single  pixel  should  be  introduced  first. 
The  figure  of  merit  for  a  charge  couple  device  (CCD)  array  should 
take  into  consideration  the  filling  efficiency  of  the  array,  spatial 
resolution  limit,  modulation  transfer  function,  quantum  efficiency 
and  integration  time  [29].  A  performance  measure  for  a  single  pixel 
may  be  used  as  a  normalization  factor  for  the  array  modulation 
transfer  function  (MTF)  response  at  zero  spatial  frequency.  As  a 
starting  point  the  voltage  responsivity  of  a  single  pixel  can  be 
expressed  as 
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[V  photon  ^  sec  cm^] 


(5-13) 


where  14  is  the  voltage  output  of  the  device,  and  Ep  is  the  uniform 
photon  irradiance  on  the  device  [photon  sec'^  cm'^]  .  If  monochromatic 
irradiance  is  considered,  the  expression  for  the  output  voltage  of 
a  pixel  can  be  rewritten  as  [29] 

=  ^p(^)  Jx  2ly  T]  e  (CE)  [V]  (5-14) 

where  AxAy  is  the  resolution  element  area,  is  the  integration 
time,  1),  is  the  quantum  efficiency,  e  is  the  electronic  charge  and 
CE  is  the  conversion  efficiency  of  the  output  preamplifier 
(V/electron) .  The  quantum  efficiency  is  defined  as  the  probability 
that  a  photoelectron  is  produced  when  a  photon  is  incident  on  the 
detector.  Substituting  in  the  responsivity  expression  Eq. (5-13) , 
the  spectral  array  photon  responsivity  can  be  expressed  as 

~E  U)  ^  ^integ  ^  ®  [V photon'^cm^sec]  .  (5-15) 

The  responsivity,  in  Eq.  (5-15)  is  valid  for  a  pixel  or 
for  an  array  at  zero  spatial  frequency.  To  include  the  resolution 
effects  would  require  that  the  responsivity  would  be  multiplied  by 
the  MTF  of  the  array  in  two  orthogonal  directions,  thus  getting  two 
expressions  for  spectral  array  photon  responsivity  in  the  x 
and  y  directions . 

The  noise  equivalent  quantum  flux  (NEQF)  is  the  counterpart  of 
the  NEP  for  a  focal  plane  array.  The  NEQF  is  defined  as 
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NEQF  . 


noise 


[photons  sec'^  cm'^] 


(5-16) 


SI 


array 

P 


where  NEQF  has  units  of  photon  sec'^  cm'^  and  noise  is  in  V. 

The  voltage  noise  of  the  device  is  caused  by  all  the 
contributors  shown  in  Fig.  5-1.  A  detailed  discussion  of  the 
detector  noise  can  be  found  in  references [22 ]  [27]  [29]  .  Hence 

noise  .  ^jlP^CE  [V]  (5-17) 

where  An„^  is  the  sum  of  all  the  noises  in  quadrature  in  units  of 
electrons  (statistical  independence) . 

Substituting  Eg. (5-15)  and  Eq.{5-17)  in  Eq.(5-16)  one  finds 

\fKP  o  1 

NEQF.  - - -  [photons  cm  sec  ]  ^'^-18'^ 

Ax  Ay  t.  ,  n 

intecr  ' 


This  is  the  photon  irradiance  flux  on  the  array  that  would  give  a 
signal-to-noise  ratio  of  1,  and  is  a  measure  of  the  minimum  photon 
irradiance  level  that  can  be  detected.  The  threshold  signal-to-noise 
level  will  still  be  required  for  detection.  The  value  of  p  will 
depend  on  the  operating  conditions. 


81 


Focal-plane  array 


state) 

•  Thermal 
generation 
(dark-current 
type) 

•  Clock 
feedthrough 


Figure  5-1.  Noise  sources  in  a  IR  seeker  using  a  CCD. 


In  the  case  of  background-limited  infrared  photodetection 
operation  (BLIP)  the  dominant  noise  is  shot  or  G-R  noise  associated 
with  the  photon  irradiance.  One  can  then  write  the  NEQF  as. 


NEQF  . 


j2G^r]E  Ax  Ay  t.  ^ 
Ax  Ay  t  .  ^  1] 

.  xnteg  ' 


(5-19) 


nEQF  -  G 


2E 


Ax  Ay  t.  ^  n 

integ  ‘ 


(5-20) 


where  G  is  the  photoconductive  gain  (G=l  for  photovoltaic  operation, 
and  0<G<1  for  photoconductor  operation)  ,  and  is  the  integration 
time . 

A  photoconductor  which  exhibits  performance  at  theoretical 
limits  is  said  to  operate  in  the  BLIP  mode.  In  this  case  the 
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tja.ckg'JTOu.nd  photon  flux  and.  the  atmosphonic  transmission  dstormin© 
the  detector  noise.  In  this  circumstance  the  flucuations  in  incident 
photon  flux  can  be  the  dominant  source  of  noise.  For  the  case  of  an 
IR  seeker  where  most  of  the  field  of  view  contains  background  and 
not  signal  photons,  and  where  the  background  photon  flux  due  to  dome 
emission  is  significantly  high,  the  term  BLIP  is  truly  appropriate. 
This  is  further  enhanced  by  the  fact  that  the  IR  sensor,  in  current 
systems,  is  cooled  well  below  ambient  temperature. 

B.  SIGNAL-TO-NOISE  RATIO  DEGRADATION  DUE  TO  DOME  EMISSION 

In  a  high  speed  flight  situation  (AL>1.2)  the  missile  heats  up 
due  to  aerodynamic  heating.  As  was  discussed  in  Chapter  I,  the 
temperature  rise  is  particularly  severe  in  the  stagnation  regions  as 
happens  with  the  IR  dome.  As  the  dome  heats  up,  it  will  emit  gray- 
body  radiation  over  the  entire  electromagnetic  spectrum.  This  will 
degrade  the  IR  sensor  performance  in  terms  of  detection  and  tracking 
ability.  The  effect  of  background- induced  noise  from  the  hot  dome  on 
the  sensor  performance  has  been  characterized  and  explained  by  some 
authors,  in  particular  Dr.  Claude  Klein.  This  section  and  the 
nonimaging  range  equation  subsection  summarizes  some  of  his  work. 

With  the  assumption  that  the  seeker  system  has  a  wavelength- 
independent  optical  efficiency  over  the  spectral  bandpass,  that  the 
sensor  subsystem  is  detector-noise  limited,  and  that  all  detector 
elements  operate  in  a  BLIP  mode,  the  peak-signal  to  rms-noise  ratio 
obtained  at  the  output  terminals  of  the  signal -processing  unit  can 
be  expressed  as  follows [4] : 
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(5-21) 


SNR  - 


H 


eff 


NEI 


Here  the  symbol  [W/cm^]  designates  the  net  effective  irradiance 

originating  from  the  target, 

H,,.  I  (X/X^)H^x^dX  (5.22) 

where  Zp  is  the  path  transmittance  and  is  the  spectral  irradiance 
originating  from  the  target  and  measured  at  the  collecting  aperture. 
NEI  is  the  noise-equivalent  irradiance. 


NEI  .  NEI^^j^yW^ 


(5-23) 


referred  to  the  peak-response  wavelength  of  the  detector  package,  Ap, 
and  expressed  in  terms  of  the  laboratory-measured  NEI^  and  the 
background  photon- flux  ratio  Under  laboratory  conditions  there 
is  no  sunlight,  and  the  system  is  in  thermal  equilibrium  (the  seeker 
and  its  surroundings  are  at  the  same  temperature) .  Under  background- 
limited  conditions  (BLIP)  ,  the  noise  voltage  is  seen  to  be 
proportional  to  the  square  root  of  the  flux  of  incident  background 
photons (Poisson  photon  statistics),  so  that  the  signal-to-noise 
ratio  (SNR)  in  the  presence  of  hot  IR  dome  radiation  can  be  derived 
from  the  "dark-system"  SNR  simply  in  the  form 


SNR  . 


SNR 
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(5-24) 
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The  ratio  thus  determines  the  performance  degradation  induced 

by  environmental  factors  and  must  be  evaluated  with  some  care.  Let 
one  consider  that  the  noise-generating  photons  originate  from  three 
distinct  sources:  the  detector  surroundings,  the  background  scene, 
and  the  hot  dome.  In  terms  of  spectral  quantum  irradiances,  the 
total  effective  photon-flux  density  O  thus  amounts  to 

*  ■  /  •  \,D  •  (5-25) 

where  the  terms  in  the  integrand  refer  to  the  spectral  photon  flux 
contributions  from  the  background,  dome  and  detector  surroundings 
respectively.  The  temperature  distribution  on  the  dome  in  practical 
designs  does  not  vary  significantly.  So  the  dome  can  be  considered 
to  be  at  a  uniform,  averaged  temperature  .  The  spectral  region  of 
interest  usually  is  narrow  enough,  implying  that  the  emittance  does 
not  change  significantly,  so  the  photon  flux  can  be  expressed  in  the 
following  compact  manner [4], 

$  (1  -  S)  Q(TJ  .  (1/T^  -  1)  Q(TJ  .  f  Q(T^) 

(J,  ■  Q  (T  )  ■  (5‘26) 

Here  7  represents  average  dome  emissivity,  is  the  transmission  of 
the  optics (about  50%  for  state-of-the-art  seekers),  Q(T)  represents 
the  effective  quantum  exitance  at  absolute  temperature  T  defined  as 

Q-  J  dX  .  [photons/ cm^-sec]  .  (5-27) 
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Under  laboratory- type  conditions,  the  effective  background  photon- 
flux  density  received  by  a  radiation-shielded  detector  amounts 
to[26] 

-K,  -  (NA)^  f  (TJ  dX  ,  (5.28) 

Here  NA  is  the  numerical  aperture  of  the  optical  system  and  tf 
represents  the  (cold)  filter  transmittance, assuming  ideal  spatial 
shielding  and  spectral  filtering. 

C.  TARGET  ACQUISITION  RANGE  DEGRADATION  DUE  TO  DOME  EMISSION 
1 .  Range  Equation  for  Non-Imaging  Systems 

A  point  target  of  spectral  radiant  intensity  gives  rise  to 
a  spectral  irradiance 


H- 


X  ■ 


X  p 


(5-29) 


at  the  collecting  aperture  of  the  seeker  system,  where  ip  represents 
the  overall  path  transmittance,  including  the  IR  dome  losses,  and  R 
is  the  target-collector  range.  Hence,  and  in  accordance  with  Eqs.(5- 
20)  and  (5-21) ,  we  have 


.  f  (X/X^)  j^T^dX 


(5-30) 


where  J^ff  refers  to  a  net  effective  apparent  target  radiant 
intensity  defined  as  follows: 
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SNR  . 


(5-31) 


NEI 

The  effective  radiant  intensity  varies  with  slant  range  and  IR  dome 
temperature  through  the  path  transmittance  Zp.  Target  acquisition  at 
a  given  probability  of  detection  and  a  fixed  false-alarm  rate 
requires  a  minimum  signal-to-noise  ratio  (SNR)^^^.  It  follows  that 
the  detection  or  acquisition  range  is  a  solution  of  Eq. (5-29) .  If  Rp 
now  designates  the  acquisition  range  for  the  "cold"  system  obtained 
under  laboratory- type  conditions,  the  degradation  induced  by  heating 
up  of  the  IR  dome  can  be  characterized  by [25] 

1/2 

'  (5-32) 


—  ^  I 

R^  NEI/ (NEI) 


(neglecting  path  transmittance  variation)  .  Using  the  relationship 
from  Eq. (5-31) 

(SNR)  (NEI) 


1/2 

(5-33) 


so  that  using  Eq. (5-21)  , 

R  ^  1 

Rp  (4>/<D^)"^^ 


(5-34) 


since  is  a  relatively  weak  function  of  both  R  and  ,  for 

scenarios  of  practical  interest. 
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FLUX  DENSITY  (3to5  photons/cm^ ■  sec  oi  deiecior) 


Figure  5-2  Photosaturation  characteristics  of  a  lead  telluride  detector 
subjected  to  background  illumination  originated  from  a  hot  magnesium 
fluoride  window [26]. 


The  dark-system  NEI  relates  to  design  parameters  in  the 
following  manner [25] : 


(NEI) ^ 


[(Af)^ 


(5-35) 


where  D*  (A^)  refers  to  the  peak  specific  detectivity  of  the  sensor 
package,  under  conditions  of  thermal  equilibrium,  at  room 
temperature,  and  is  an  electric  response  factor. 

Figure  5-2  shows  calculated  and  computed  minimum  detectable 
signal  power  relative  to  minimum  detectable  signal  power  in  the 
absence  of  window  radiation [26 ] .  This  example  shows  the  strong 
effect  of  the  dome  emission  on  detection  and  tracking. 
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2 .  Range  Equation  for  Imaging  Systems 

The  NETD  Eq. (5-10)  is  valid  either  for  a  BLIP  or  a  non-BLIP 
detector.  In  this  section  the  NETD  for  a  BLIP  detector  irradiated  by 
a  hot  window  will  be  derived.  As  will  be  shown  the  range  equation 
for  an  imaging  system  will  be  NETD  dependent. 

The  theoretical  limit  for  D*  for  a  photoconductor  is 


D-  (X) 


X 

2hc 

\ 


1/2 


(5-36) 


where  rj^  is  the  quantum  efficiency  and  is  the  background  quantum 
flux  density  [photons /cm^-sec]  incident  on  the  detector.  The  severity 
of  the  IR  seeker  signal-to-noise-ratio  limited  application  dictates 
that  a  cold  filter,  with  appropriate  spectral  response,  should  be 
placed  immediately  preceding  the  detector.  The  is  given  then  by 


/  (T,)  dX  , 


(5-37) 


where  it  assumed  that  the  detector  is  a  Lambertian  receiver  with  an 
effective  collecting  solid  angle  of  n,  and  where  Q^s  is  the  ideal 
effective  angle  to  which  the  detector  cold  shield  reduces  reception 
of  background  radiation,  and  t^o  is  the  transmittance  of  the  cold 
filter.  At  this  point,  it  is  convenient  to  define  the  several 
individual  transmittances  that  affect  the  sensor  performance.  The 
total  spectral  optical  transmittance  can  be  decomposed  in  the 
following  product 
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''■o,X  ■  ^CO,X  '  ^  UC,X  ’ 


(5-38) 


where  Tjjc,a  is  the  spectral  transmittance  of  the  uncooled  optics  and 
'^d.aC^d)  is  the  average  spectral  transmittance  of  the  dome  at  an 
averaged  absolute  temperature 

In  this  analysis  it  is  assumed  that  the  detector  is  a 
photoconductor  having  theoretical  performance  such  that 


Ad-  (A  ) 

D"  (A)  .  - - — —  for  A  £  A 

X  ^  (5-39) 

-  0  for  A  >  A  . 


Substituting  for  in  the  D*  Eq.  (5-3  6)  gives 


D-  (K) 


A 

2hc 


n 


i/i 


ri 


1/2 

g 


/  (T,)  dX 


(5-40) 


A  normalized  detectivity  parameter,  D*\j^^p(A)  ,  for  unity  quantvim 
efficiency  and  for  an  effective  cold  shielded  angle  of  n  steradians 
is  defined  by [22] 


^BLIP  • 


A 

2hc 


1 


f  ^CO  2x  (Tp)  dX 


(5-41) 


Expressing  D*  (A)  in  terms  of  D**  (A) 
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D-(K)  . 


n 


1/2 


D'slip 


rx; 


(5-42) 


Substituting  for  D*  (Ap)  in  the  NETD  equation  yields 


NETD 


(5-43) 


Here  AM/ AT  is  defined  as 


Am 

At 


/ 


D-(X) 


ST  D- (X^) 


dX 


(5-44) 


where 


d  T 


c_  /XT 


■ 


\  m2  ,  ^2 

X  (e  ^ 


1) 


(5-45) 


gives  the  effective  change  in  spectral  radiant  emittance  with 
temperature.  AT  is  the  temperature  difference  between  the  target  and 
the  background.  Using  the  relation  in  Eq.(5-39),  Eq. (5-44)  can  be 
simplified  to 


At  '  I  dT  X 

K  p 


(5-46) 


The  cold  shield  efficiency  77^5  is  defined  in  terms  of  the  actual  cold 
shield  angle  Op,  as 

(5-47) 
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where  Q^s  =  a-b  /  f^. 

The  perfect  cold  shield  is  defined  as  one  which  limits 
reception  of  background  radiation  to  the  cone  defined  by  the  optical 
system  aperture  ratio,  F„  (f /number)  .  Thus  for  an  optical  system  with 
circular  aperture  of  diameter  and  focal  length  f 


Q  . 


p 


(5-48) 


and  substituting  in  Eq. (5-43)gives[22] 


NETD  . 


_ ^  ^/^ 


Z\Af 

AT 


(5-49) 


Applying  the  same  definition  to  D*  gives 

D\  (X)  .  2  Fn  rrics  .  (5.50) 

The  Eq.  (5-49)  can  be  rewritten 


NETD  - 


D  X  D-  (X  )  — 

y  ^  o  o  *  p' 


(5-51) 


To  account  for  the  effect  of  the  hot  dome  on  the  NETD  the 
background  quantum  flux  density  from  the  dome  emission  should  be 
added  to  the  background  quantum  flux  density  from  the  scene 
background  in  the  D*  Eq.  (5-36)  .  With  the  cold  shield,  the  background 
quantum  flux  density  that  irradiates  the  detector  has  two 
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contributions : 


(5-52) 


where  is  the  quantum  flux  density  originated  by  a  uniform 
blackbody  background  corresponding  to  the  scene  at  temperature  Tgg, 
Qg  is  the  portion  of  the  background  quantum  flux  density  originated 
by  the  dome  with  a  uniform  mean  temperature  Tg  .  These  quantities  are 
defined  as  follows  [25] : 


Q 


BB 


,  /Xt_, 


(5-53) 


A^ 


c-  /XT* 


1) 


dA 


(5-54) 


with  radiometric  constant  C2  as  defined  before  and 
1 . 8837xlO^VniVsec-cm2,  e  is  the  mean  dome  emissivity  with  an 

averaged  temperature,  and  y\.2  are  the  cut-on  and  cut-off 

wavelengths  of  the  filter.  Substituting  in  Eq. (5-38)  yields 


D-  (X). 


A 

2nc 


n 


\l/2 


n 


1/2 


,X 


;  dX 


1/2 


(  Q 


BBrK 


(5-55) 


which  can  be  used  to  express  NETD  from  Eq.  (5-50)  in  terms  of  photon 
flux  components. 

Typically  the  optical  system  installed  in  the  seeker  is  a 
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modified  Cassegrain  telescope.  This  will  introduce  some  energy 
losses  characterized  by  the  circular  aperture  attenuation 
coefficient,  K(cir) ,  and  a  Cassegrain  obscuration  coefficient 
K(cas) .  If  these  losses  are  considered,  the  general  NETD  becomes 


NETD  =  ^ 


1/2  Fn  (Qcs)^^^ 


yfaP  K(cas)K(cir)  [  n. 


1/2 


J  ^  BB  A  ^BB  ^  ‘^^D,  X  ( 


1/2 


f 


7  ^^2  /^I’bB 

M  ^ _ 

2  nBfA.  p,  /3rp 

Tbb  e  -  1 


dl 


(5-56) 


where  x^c  is  the  cold  optics  transmittance  (under  laboratory 
conditions),  x^  is  the  optics  transmittance  (during  operation),  M^gi 
and  Mp  ,,  are  the  spectral  radiant  exitances  of  the  background  and  of 
the  dome  respectively.  Expanding  the  terms  of  the  integrands  of 
Eg.  (5-56)  another  format  to  the  general  NETD  equation  can  be  found. 


NETD 


yfaPD^v^  K(cas)K(cir)  (  »7,j  j 


/  i  p 


]  dX 


1/2 


f- 


3  ^2  /I'^BB 


3  C_2  /XTqq 


dX 


(5-57) 

The  last  two  expressions  for  the  NETD  take  into  consideration  the 
influence  of  radiation  by  an  aerodynamic ally  heated  dome  on  the 
sensor.  The  added  "background" -induced  noise  will  lower  the 
effective  performance  of  the  system  in  terms  of  detection  range  and 
imaging  capability. 

The  signal-to-noise  ratio  at  maximiom  detection  range  is 
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(SNR) 


AT  -T  (R) 

_ p  ' 

NETD 


(5-58) 


where  Tp(R)  is  the  total  atmospheric  path  transmission  for  the  system 
spectral  bandpass  and  R  is  the  detection  range.  Here  AT  is  the 
differential  temperature  of  the  target  with  respect  to  background. 
If  AT  is  normalized  with  respect  to  the  SNR,  defining[30] 

^  ^  Differential  Temperature  _  AT  (5-59) 

^  (SNR)  ~  (SNR) 

a eg  '  aeg 

and  substituted  in  Eq. (5-58)  then  after  rearranging  one  obtains 


T^(R)  . 


NETD 


(5-60) 


To  solve  for  range,  the  total  transmittance  of  the  optical  path 
should  be  determined  first.  For  a  given  set  of  meteorological  and 
oceanographic  conditions,  the  transmittance  falls  exponentially  with 
range,  as  follows 


x^(R)  .  ae  (5.61) 

where  a  and  (3  are  constants  to  be  determined.  Substituting  Eq.  (5- 
61)  in  Eq. (5-60) ,  gives 
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D.  DETECTOR  NONUNIFORMITY  NOISE 

This  noise  is  associated  with  the  variation  in  output  from 
pixel  to  pixel,  or  the  variations  across  the  entire  array.  It  can  be 
caused  by  two  effects: (a)  fixed  pattern  noise  and  (b)  responsivity 
variations.  The  influence  of  nonuni foirmity  on  FPA  performance  may 
impair  the  detector  of  achieving  theoretical  (BLIP)  sensitivity. 

Nonuniformity  limitations  on  the  NETD  in  many  systems  are  more 
damaging  than  temporal  noise  since  a  pattern-noise-limited  image 
cannot  be  improved  by  frame-to- frame  averaging,  whereas  the 
sensitivity  of  a  temporal-noise-limited  imager  can  be  improved  by 
such  averaging. 

1 .  Fixed  Pattern  Noise 

Fixed  pattern  noise  is  a  nonuniformity  appearing  across  the 
array  which  does  not  change  with  irradiance [29 ] .  Several  effects 
contribute  to  fixed  pattern  noise.  Themtially  generated  carriers  are 
not  uniform  across  the  array.  Local  defects  in  the  semiconductor 
crystal  lattice  cause  a  variation  in  dark-current  generation  across 
the  array. 

A  major  contributor  to  fixed  pattern  noise  is  the  errors  in  the 
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fabrication  process .  The  errors  introduced  by  inaccurate  pixel  areas 
caused  by  masking  misalignments  can  lead  to  10%  nonuniformities.  The 
fixed  pattern  noise  is  usually  compensated  by  image  processing. 

2  .  Responsivity  Variation 

Responsivity  variation  is  a  nonuniformity  which  has  a  variable 
effect  as  a  function  of  irradiance .  The  responsivity  variation  from 
pixel  to  pixel  is  determined  by  doping  variations  across  the  device. 
Compensation  for  responsivity  is  .more  difficult  then  for  fixed 
pattern  noise-  requiring  "two  point"  calibration  as  a  minimiam. 

3  .  Residual  Nonuni foinnity 

After  compensation  the  residual  non-uniformity  can  be 
considered  as  a  source  of  noise.  This  noise  can  be  considered  as  an 
added  percentage  of  photon  noise,  because  it  is  considered  to  be 
correlated,  in  the  following  manner, 

^  (Qbb  •  (5-64) 

where  n  is  the  rms  value  of  the  noise  in  volt,  g  is  the  non¬ 
uniformity  coefficient,  Q^b  and  are  the  background  photon  fluxes 
on  the  detector  from  the  scene  and  seeker  dome  as  defined  before  by 
Eq.  (5-53)  and  Eq.  (5-54)  respectively.  Array  non-uniformity  is  a 
problem  area  that  has  been  discussed  extensively  in  the  literaature. 
For  state-of-the  art  sensors  the  array  non-uniformity  coefficient 
varies  between  0.01  and  0.00 01 [30]. 

E.  COMPUTATION  OF  THE  DETECTION  RANGE 

Acquisition  range  for  an  IR  seeker  of  an  air-to-surface  anti¬ 
ship  supersonic  missile  with  a  256x256  CCD  staring  sensor  was 
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E.  COMPUTATION  OF  THE  DETECTION  RANGE 

Acquisition  range  for  an  IR  seeker  of  an  air-to-surf ace  anti¬ 
ship  supersonic  missile  with  a  256x256  CCD  staring  sensor  was 
computed  for  a  midlatitude  scenario  (LOWTRAN  6)  ,  in  winter  and 
summer,  for  two  background-to-target  normalized  differential 
temperatures  (lCi,=2 . 5  and  Xj,=4.0)  and  for  two  missile  flight  profiles, 
level  flight  at  20  m  and  3km  altitude,  and  for  a  range  of  mean  dome 
temperatures  of  300,  600,  900,  1043  and  1200K,  corresponding  to 

missile  flight  speeds  up  to  Mach  4. 

The  detection  or  acquisition  range  was  computed  using  Eq. (5-63) . 
To  compute  the  NETD  a  MS-FORTRAN  program  developed,  by  S.T.  Smith 
(Naval-Air  Warfare  Center,  China  Lake)  was  used.  The  program 
computes  the  NETD  for  different  detector  materials,  integration 
times,  dome  materials,  dome  temperature,  array  nonuniformities,  fill 
efficiency,  background  temperature,  collector  aperture  and  optical 
configurations .  The  NETD  was  computed  for  each  set  of  parameters 
shown  in  Table  5-1,  with  LOWTRAN  atmospheric  transmission. 

The  total  transmittance  of  the  optical  path  was  determined 
first.  The  LOWTRAN6  computer  code  was  used  to  predict  atmospheric 
transmittance.  Figure  5-3  shows  the  total  computed  transmittance 
for  a  slant  path  from  3km  altitude  to  sea  level  and  for  a  horizontal 
path  at  20  meters  altitude  for  several  ranges  in  a  Navy  Maritime 
environment (see  Appendix  B)  in  the  3. 5-4. 2  ym  and  8.0-10.5  ym,  for 
the  midlatitude  summer  and  winter  models. 
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Figure  5-3.  LOWTRAN  6  computed  atmospheric  transmittance.  *  3.5-4.2pm  +  8.0-10. 5pm. 
Navy  maritime  aerosol  model.  (a)(b)  Midlatitude  Summer ;visibility  44.2  km  ;Humidity  76.2%; 
(c)(d)  Midlatitude  Winter;visibility  43.9km; Humidity  77.1%). 

Using  a  least-squares  fit,  the  coefficients  a  and  3  of  Eq. (5- 
61)  were  computed  for  the  four  different  scenarios  for  the  two 
spectral  bands.  The  curve-fit  parameters  for  atmospheric 

transmission  and  the  computed  NETD  for  the  each  set  of  sensor 
parameters  (see  Table  5-1)  and  dome  temperatures  of  300,  600,  900, 
1043  and  1200K,  were  used  to  compute  the  seeker  detection  range  by 
applying  Eq. (5-63) .  The  selected  dome  materials  were  the  sapphire 
and  the  diamond  for  the  3.5-4.2iam  band  and  8.0-10.5iam  spectral  band 
respectively.  These  materials  were  matched  to  the  spectral  bands 
where  they  exhibit  better  transmission  properties.  These  materials 
are  considered  to  withstand  the  severe  thermal  environments  of 
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advanced,  highly  supersonic  missiles.  The  temperature  of  1144K  was 
found  in  Chapter  IV  to  be  the  mean  dome  temperature  of  the 
conceptual  missile  at  Mach  4  at  sea  level. 

The  normalized  differential  temperatures,  K^2.5  and  K^A.O  were 
chosen  for  a  SNR  of  2.0.  This  implies  that  we  are  looking  at 
temperature  differences  of  5 . OK  and  8. OK  between  the  target  and  the 
background . 

F.  RESULTS  AND  DISCUSSION 

Figure  5-4  through  Fig. 5-7  show  the  computed  NETD  for 
"residual"  array  nonuniformity  coefficient  after  compensation  of 
0.1,  0.01,  0.001  and  0.0001  respectively  for  the  two  selected  dome 
materials.  Fig. 5-8  through  Fig. 5-23  show  the  computed  detection 
ranges  for  the  seeker. 

Before  starting  a  detailed  analysis  of  the  detection  range 
computation  results  it  should  be  noticed  that  the  values  of  the 
computed  NETD  depend  on  the  area  of  the  collector  aperture (see 
Eq.(5-10)).  So  all  the  discussion  relates  to  the  nose  shape  and 
dimensions  of  the  hemispherical  nose  configuration  of  Chapter  IV. 
The  quantitative  results  obtained  cannot  be  immediately 
extrapolated  to  other  seekers .  However  the  results  obtained  show  the 
trends  and  influence  of  the  several  system  parameters  on  the  IR 
seeker  performance.  A  LRASM  (Long  Range  Anti-Ship  Missile)  would 
enable  the  installation  of  a  larger  primary  optical  aperture 
resulting  in  better  detection  range  and  good  resolution  for  imaging. 

Analyzing  the  results  it  is  clear  that  the  array  nonuniformity 
has  strong  effect  on  the  seeker  performance.  The  array  nonuniformity 


100 


is  particularly  important  at  high  supersonic  speeds  when  the  dome 
temperatures  become  particularly  high.  For  the  case  of  a  Mach  4 
missile  the  seeker  becomes  blind  (detection  range  falls  to  zero)  for 
residual  array  nonuni fointii ties  of  10  percent.  For  the  other  cases  of 
improved  array  nonuniformity  the  detection  range  of  the  system 
improves  considerably.  The  sapphire  dome  systematically  shows  a 
better  performance  than  the  diamond  dome  in  terms  of  detection 
range . 

In  terms  of  environmental  influence  on  the  performance  the 
winter  condition  with  a  slant  path  scenario  gives  the  best  results 
in  opposition  to  the  summer  scenario  for  a  sea-skimmer  missile  at 
an  altitude  of  20m. 

The  results  show  clearly  that  the  designer  should  determine 
carefully  the  effect  of  array  nonuniformity  and  the  dome  temperature 
reached  during  flight  by  the  seeker. 

The  model  used  does  not  take  into  account  target-to-background 
contrast  reflectance.  A  better  model  should  include  this  effect  if 
a  better  target-to-background  discrimination  is  desired.  This  effect 
is  expected  to  be  larger  in  the  MWIR  spectral  band  than  the  LWIR 
spectral  band. 

For  aim-point  considerations  a  minimum  of  100  pixels  on  target 
is  desired  (Johnson  criteria) .  To  calculate  the  aim-point  tracking 
range  one  should  take  into  account  the  geometrical  characteristics 
of  the  target  and  the  instantaneous  field  of  view  of  the  seeker 
optical  system. 
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FPA. 


256x256,  30  |jm  pitch 


Aperture  Diameter . 2.25  cm 

FOV . 4.46° 

IFOV . 0.2° 

Diffraction  limit  blur . 0.13mrad 

Obscuration  coefficient . 0.1 

Dome  Material 

3.5  -  4.2pm  band  (  Sapphire  ) 

8.0  -  10.5pm  band  ( ZnS  and  Diamond  ) 

Detector  Material 

3.5  -4.2pm  band  . InSb 

8.0  -  10.5pm  band . HgCdTe 

Array  Nonuniformity 

0.001,  0.01, 0.10,1.00 

Dwell  Time 

Staring . 4.0  ms 

Dome  shapes 

hemispherical 

elipsoidal 


Table  5-1.  Seeker  Parameters.  Study  Cases. 
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NETD,  [K]  NETD,  [K] 
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100 


101 


256x256  FPA  Integ.  Time:  4iiisec 
1  percent  Nonunif. 


Diamond 
8.0-10.5  microns,.. 


^/Sapphire 
3.5-4.2  microns 


400  600  800  1000 
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Figure  5-4 
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Figure  5-6 
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Figure  5-7 


Detection  Range,  [km]  Detection  Range,  [km] 


256x256  FPA  Integ.  Time:  4  ms  Kt=4.0 


256x256  FPA  Integ.  Time:  4  ms  Kt=4.0 
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256x256  FPA  Integ.  Time:  4  ms  Kt=  4.0 


256x256  FPA  Integ.  Time:  4  ms  Kt=  4.0 
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Detection  Range,  [km]  Detection  Range,  [km] 


256x256  FPA  Integ.  Time:  4  ms  Kt-4.0 


256x256  FPA  Integ-  Time:  4  ms  Kt=4.0 
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Detection  Range,  |kin]  Detection  Range,  [km] 


256x256  FPA  Integ.  Time:  4  ms  Kt=  4.0 


256x256  FPA  Integ.  Time:  4  ms  Kt=4.0 
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Figure  5-15 
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Detection  Range,  [km]  Detection  Range,  [km] 


256x256  FPA  Integ.  Time:  4  ms  Kt=  4.0 
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256x256  FPA  Integ.  Time:  4  ms  Kt=4.0 


256x256  FPA  Integ.  Time:  4  ms  Kt=4.0 
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Detection  Range,  [km]  Detection  Range,  [km] 


256x256  FPA  Integ.  Time:  4  ms  Kt=  4.0 


256x256  FPA  Integ.  Time:  4  ms  Kt=4.0 
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Dome:  Sapphire  (  e^lxliX^ ) 

Detector:  InSb 

Spectral  Band:  3. 5 -4. 2  pm 


Dome:  Diamond  (  e^SxlQr^ ) 

Detector:  MCT 

Spectral  Band:  8. 0-10,5 pm 


Dome 

Temperature 

[K] 

NEI 

[W/cm^] 

NETD 

[K] 

NEI 

[W/cnP] 

NETD 

[K] 

300 

3.293x10-'^ 

2.236x10-^ 

3.545x10-'" 

1.026x10-' 

600 

3.958x10'*^ 

2.689x10-' 

5.606x10-'" 

1.623x10' 

900 

7.335x10'^ 

4.982x10-' 

9.336x10-'" 

2.702x10-' 

1043 

9.705x10*^ 

6.592x10-' 

1.140x10-" 

3.301x10' 

1200 

1.270x10-'" 

8.628x10-' 

1.379x10-" 

3.993x10-' 

Table  5-2.  Computed NEI and NETD  for  a  256x256  FPA  with  0.01%  Nonuniformity. 


Dome:  Sapphire  (  €=lxlP^) 
Detector:  InSb 

Spectral  Band:  3. 5-4. 2  pm 

Dome:  Diamond  (  e~5xl(T^ ) 
Detector:  MCT 

Spectral  Band:  8.0-10.5 pm 

Dome 

NEI 

NETD 

NEI 

NETD 

Temperature 

[K] 

[W/cm^J 

[K] 

[W/cm’J 

[KJ 

300 

6.4471x10-'^ 

4.3790x1(1^ 

3.2358x1(1” 

9.3677x10’ 

600 

9.0596x1 

6.1534x1(1’ 

5.2919xia” 

1.5320x10’ 

900 

2.7617x1(1’^ 

1.8758x10-’ 

9.0180x10-” 

2.6107x10’ 

1043 

4.4396x1 Or’^ 

3.0155x10-’ 

1.1082x10” 

3.2083x10’ 

1200 

6.8278x1(1” 

4.6375x1 a’ 

1.3473x10” 

3.9003x10’ 

Table  5-3.  Computed  NEI  and  NETD  for  a  256x256  FPA  with  0.1  %  Nonuniformity. 
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Dome:  Sapphire  ( e=lxl(T^) 

Detector:  InSb 

Spectral  Band:  3. 5 -4. 2  pm 


Dome:  Diamond  (  e=5xlB^ ) 

Detector:  MCT 

Spectral  Band:  8.0-10.5 pm 


Dome 

Temperature 

[K] 

NEI 

fW/cm^] 

NETD 

[K] 

NEI 

[W/cm^J 

NETD 

[K] 

200 

3.789x10-' 

9.358x10-' 

600 

5.570x10-' 

1.531 

900 

1.818 

2.609 

1043 

2.9574 

3.207 

1200 

4.5795 

3.870 

Table  5-4.  Computed  NETD  for  a  256x256  FPA  with  1.00  %  Nonuniformity. 


Dome:  Sapphire  (  e=lxl0'^) 
Detector:  InSb 

Spectral  Band:  3. 5 -4. 2  pm 

Dome:  Diamond  (  e=5xl(T^ ) 
Detector:  MCT 

Spectral  Band:  8.0-1 0.5 pm 

Dome 

NEI 

NETD 

NEI 

NETD 

Temperature 

[K] 

[W/cmf 

[K] 

[W/cmf 

[K] 

300 

3.2933x10-'^ 

2.2368x10-' 

3.5448x10-''' 

1.0262x10-' 

600 

3.958exl0-'’ 

2.6886x10-' 

5.6063x10-''' 

1.6230x10-' 

900 

7.3345x10-'^ 

4.9817x10-' 

9.3362x10-''' 

2.7028x10-' 

1043 

9.7048x10-'^ 

6.517x10-' 

1.1401x10-" 

3.3007x10-' 

1200 

1.2703x10-''* 

8.6283x10' 

1.3779x10-" 

3.9930x10-' 

Table  5-5.  Computed  NEI  and  NETD  for  a  256x256  FPA  with  0.0001  %  Nonuniformity. 
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Dome:  Sapphire  (  €=lxl(T^) 
Spectral  Band:  3. 5 -4. 2  pm 

Dome:  Diamond  (  e=5xl0'^) 
Spectral  Band:  8. 0-1 0.5 pm 

Dome 

Computed  Photon  Incidence 

Computed  Photon  Incidence 

Temperature 

[K1 

[Photons/sec  cm^  sr] 

[Photons/sec  cm^  sr] 

300 

7.9172x10'" 

1.1486x10" 

600 

1.6141x10'^ 

1.8792x10'^ 

900 

3.8033x10" 

3.2031x10" 

1043 

6.1885x10'^ 

3.9366x10" 

1200 

9.5829x10" 

4.7860x10" 

Table  5-6.  Computed  Photon  Incidence  in  the  focal  plane. 
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VI.  COMPUTATION  OF  PHOTON  NOISE  FROM  HOT  DOME 

The  characterization  of  fixed  pattern-noise  is  an  important 
aspect  in  evaluating  the  performance  of  a  staring  array.  The  hot  IR 
dome  will  induce  on  the  sensor  a  background  photon  noise 
distribution  that  can  be  considered  like  a  fixed  pattern-noise.  In 
discussing  the  "hot  dome  problem"  the  question  arises  as  to  what 
dome  temperature  distribution  should  be  used  in  the  calculation  of 
the  dome  emission  effects  on  a  staring  sensor  performance.  Is  it 
sufficient  to  consider  the  temperature  of  the  dome  to  be  uniform, 
and  use  a  mean  temperature,  or  will  the  background  photon  flux 
induced  by  dome  emission  impose  a  significant  fixed  pattern  noise  on 
the  detector  array?  A  useful  simplification  of  the  problem  would  be 
to  consider  the  photon  flux  constant  through  the  detector  area  and 
consider  an  averaged  dome  temperature . 

In  this  chapter  these  questions  will  be  investigated  by 
computing  the  radiative  interchange  between  the  hot  dome  and  the 
detector  array. 


A.  RADIATION  TRANSFER 

Given  the  temperature  distribution  on  the  nose,  we  are  able  to 
compute  the  IR  emission  due  to  the  aerokinetic  heating  of  the 
window.  In  our  case  we  want  to  compute  the  self -emission  generated 
by  the  hot  window  that  will  fall  on  the  seeker's  IR  sensor.  As 
assumptions,  we  will  consider  that  the  window  has  reached  thermal 
equilibrium,  and  so  the  temperature  of  the  dome  is  the  recovery 
temperature,  that  the  window  is  thermally  thin,  and  that  the  dome 
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inner  surface  is  not  coated  with  any  material.  In  order  to  compute 
the  window's  self-emission  sensed  by  the  detector,  an  optical  system 
must  be  assiimed.  In  this  case  a  small  telescope,  composed  simply  of 
a  single  germanium  lens,  is  adopted.  This  optical  design  may  seem 
unrealistic  (no  field  lens  or  light  pipe  are  used),  but  in  fact,  it 
can  serve  as  a  benchmark  to  compare  some  of  the  aero-thermal  effects 
that  degrade  the  system  image  quality  due  to  dome  emission. 

To  start  the  radiant  transfer  analysis,  the  radiant  power 
reaching  an  element  of  area  dA^  of  the  detector  surface  due  to  the 
source  element  dA„  on  window  surface  should  be  computed  first. 
For  simplicity  one  should  assume,  for  the  moment,  that  no  lens  is 
present  between  the  self-emitting  window  and  the  detector.  The  line 
joining  the  element  areas,  of  length  makes  angles  of  e„  and  6^ 
with  the  respective  normals  to  the  surface  elements,  as  shown  in 
Fig.  6-1.  The  radiant  power  is  d^W^„,  a  second  order  differential 
because  both  the  source  and  receptor  are  elemental  areas [11] [21],  so 
that 


LdA  dA.cosd  cosB^ 
- 

dw  2 


(6-1) 


where  L  is  the  radiance  at  dA„. 

The  total  radiant  power  at  the  entire  second  surface  due  to  the 
entire  first  surface  is,  by  integration. 


LcosQ  cosQ.dA^dA 

w  d  d  w 


2 


(6-2) 


E 


Figure  6-1.  Configuration  for  interchange  between 
two  infinitesimal  elements. 

In  adding  power  rather  than  amplitudes  in  this  integration,  we 
have  tacitaly  assumed  that  the  radiation  source  emits  incoherent 
radiation  and  that  the  radiation  emitted  from  any  point  of  the 
surface  is  diffusely  distributed,  obeying  Lambert's  law.  If  we 
discretize  the  window  and  the  detector  in  finite  area  elements,  and 
consider  each  element  to  be  isothermal,  the  total  irradiance  into 
the  detector  area  element  originated  by  the  window  is 

j  L  .cosQ  cos0^  A  A_, 

•'.-E  ^ (6.3) 


where  the  window  is  discretized  in  J  finite  area  elements  and  the 
detector  in  K  finite  area  elements .  The  radiance  in  the  spectral 
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region  of  interest  of  the  window  finite  element  j  is 


L.  (X,T)- 


^  €.  (X,  T)  \ 

*  2nc'' 

/  \ 
I 

n  j 

J 

ch/XJcT .  ^ 

(  e  ^-1) 

dX 


(6-4) 


where  c  is  the  speed  of  light  in  vacuum,  h  is  Planck's  constant  and 
k  is  the  Boltzmann's  constant,  and  £’(X,T)  the  emissivity  of  the 
infinitesimal  element. 


1 


Figure  6-2.  Configuration  for  interchange  between  two  finite 
surfaces. 

As  was  shown  in  the  previous  chapter,  when  the  interest  is  in 
the  limits  imposed  by  dome  heating  on  seeker  performance,  a  more 
relevant  quantity  is  the  photon  exitance  Q  (photons  s'^  cm'^)  from  the 
window.  The  spectral  photon  exitance,  (photons  s'^  cm’^Vim’^  )  ,  is 
obtained  by  dividing  Planck's  law  by  hc/A, 
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Q^(Tr^) 


^  2nc' 

[  -  ^  ] 

J 

(6-5) 


The  integrated  photon  exitance,in  a  given  wavelength  band,  is  then 
obtained  by  integrating  Eq. (6-5)  between  the  cut-on  and  cut-off 
wavelengths , 


Q(T)  .  I  £3^  dX  . 


(6-6) 


Numerical  approximation  of  this  integral  over  a  narrow  wavelength 
range  yields 


Q(T) 


— — - O'.K' 

A  4  hc/K  kT  ^  ^  ^ 


(6-7) 


where  is  the  average  wavelength  within  the  band  Ai  and  Aj. 

Then  the  photon  flux  that  irradiates  a  finite  area  element  k  of 
the  detector  is 


*.■5: 


j  E.(X,T)Q  cosQ  cosQ.  a  a. 
£ _ J  ”j  “i* 


(6-8) 


where  is  the  black  body  exitance  of  the  element  j. 
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B.  COMPUTATION  METHOD 

Now  let's  consider  that  between  the  window  and  the  detector 
there  exists  a  thick  converging  lens.  This  system  is  not  in  thermal 
equilibrium.  In  fact  for  simplification  consider  that  the  detector 
and  lens  are  cooled  to  zero  absolute.  To  compute  the  irradiance  on 
the  detector  due  to  the  thermal  emission  of  the  heated  window,  a  ray 
tracing  must  be  done.  In  this  method,  the  paraxial  system  matrix  is 
used  to  compute  the  image  surface  for  the  lens  conjugate  to  the  hot 
irdome  (infrared  dome)  as  object  surface.  This  image  surface  is  then 
treated  as  an  equivalent  input  surface  viewed  directly  from  the 
detector  plane.  The  conjugate  surface  is  found  by  mapping  several 
points  of  the  window  by  means  of  a  matrix  transformation.  This 
virtual  surface  is  then  divided  into  diffusely  emitting  finite 
elements.  Radiative  transfer  calculations  between  each  detector 
element  and  the  finite  elements  of  the  image  surface  are  carried 
out,  taking  into  consideration  the  respective  shape  factors.  The 
contributions  of  all  finite  elements  are  integrated,  enabling  the 
computation  of  the  photon  flux  that  irradiates  each  photo-sensitive 
detector  element.  This  procedure  is  applied  to  all  the  elements  of 
the  detector  array.  An  array  fill  factor  of  1.0  is  assumed. 

1 .  System  Ray-Transfer  Matrix 

If  we  consider  the  temperature  field  of  the  window  to  be 
axisymmetric,  resulting  from  zero  incidence  of  the  missile,  and  the 
optical  system  to  be  looking  in  the  missile  longitudinal  axis 
direction,  the  problem  of  computing  the  photon  flux  that  falls  into 
the  detector  becomes  a  2 -dimensional  calculation  with  cylindrical 
symmetry.  Considering  the  optical  system  consisting  of  a  single 
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lens,  a  matrix  can  be  constructed  representing  the  action  of  the 
thick  lens  on  the  radiation  emitted  from  the  window,  as  follows: 


M  -  S3  at,  S3  31^  s,  (6.9) 

where  M  the  composite  transfer  matrix  is  formed  by:  translation 
matrix,  3  1  ,  representing  the  optical  path  from  the  input  plane  to 
the  front  face  of  the  lens;  refraction  matrix,  ,  representing  the 
refraction  process  at  the  lens  front  surface;  3  2  /  translation 
matrix  through  the  lens  medium;  9^2  -  describing  the  refraction  at 

the  lens  back  surface;  and  finally  3  3  ,  describing  the  final 
optical  path  from  the  lens  back  surface  to  the  output  plane, 
conjugate  to  the  input  plane. 

The  system  matrix  has  the  following  format, 


m^i 

^12 

^21 

^22 

(6-10) 


If  the  output  plane  corresponds  to  the  conjugate  to  the  input 
plane,  the  matrix  element  must  be  identically  zero  [21]  .  Then 
is  the  transverse  linear  magnification. 

2  .  Calculation  of  the  System  Ray- Transfer  Matrix 
In  this  section  the  system  matrix  will  be  calculated  by 
multiplying  consecutively  the  several  translation  and  refraction 
component  matrices  for  paraxial  rays. 
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(6-1  la) 


(6-1  lb) 


(6-1  Ic) 


(6-1  Id) 


3 


3 


1  X 
0  1 


(6- lie) 


where  is  the  distance  from  the  input  surface  to  the  lens,  and 
i?2  are  the  radii  of  curvature  of  the  front  and  back  surfaces  of  the 
lens,  n  and  n_j  the  indexes  of  refraction  of  the  air  and  of  the  lens 
respectively,  t  represents  the  lens  thickness,  and  finally,  x  is  the 
distance  from  lens  to  the  output  plane.  Carrying  out  the  matrix 
multiplication  as  in  Eq. (6-9)  yields  the  following  matrix  elements 
and  of  fh®  system  matrix  M: 
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when  m]^2  ~  0  ,  the  output  plane  is  the  image  plane  conjugate  to  the 
input  plane.  Using  Eq. (6-12a)  the  location  x  of  the  output  plane  can 
be  then  computed: 


(6-13) 

Substituting  this  value  of  x  in  Eq. (6-12b)  gives  the  value  of 
linear  magnification  at  that  position. 

In  order  to  determine  the  system  limiting  aperture  (footprint  at 
the  dome)  the  other  elements  should  be  calculated  as  well,  leading 
to  the  following  values: 
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3 .  Computation  of  the  Footprint  at  the  Dome 

The  computation  of  the  footprint  at  the  dome  of  the  system  is 
necessary  in  order  to  determine  which  portion  of  the  dome  will 
irradiate  the  detector,  and  so  contribute  to  the  background  noise. 
This  can  be  done  by  tracing  back  the  limiting  optical  rays  from  the 
detector  that  pass  through  the  lens  aperture.  Again  considering  the 
same  system  geometry,  the  problem  becomes  2 -dimensional .  We  have 


y 

yo 

sinOC 

.  M' 

sinOL 

0 

(6-15) 


where  is  the  angle  subtended  by  the  ray  from  the  detector  passing 
through  the  lens  margin,  as  shown  in  Fig.  6-4,  and  the 
composite  matrix  M'  defined  as  follows: 

M  .  ^2  ^1  (6-16) 

where  ,  32'  =  ^2  ,  9?2 '  =  5R2  >  and  the  translation  matrix 
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from  the  detector  into  the  lens  back  surface  is  defined  as  follows: 

(6-17) 

The  matrix  multiplication  as  written  in  Eq.{6-14)  leads  to  the 
following  matrix  M'  elements 


1) 


(6-1 8a) 


(6-1 8b) 


(6-1 8c) 


(6-1 8d) 


Substituting  the  following  values  into  the  respective  expressions. 


n-1.000  t=0.300  cm 

nL(lOiim)  =  4.003174  =  9.09099  cm 

Ri  =  -50.000cm  =  1.500  cm 

R2  =  60.000  cm 
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yields  the  following  matrix: 


M'  . 


0.996249 
-5.  683183e-2 


9.131830 
-1. 483700 


(6-19) 


with  given  by 


a 


o 


tan 


(6-20) 


Using  Eg. ( 6-13 )  , 


y 

0*996249 

9.131830 

0 

sinOf 

-5. 683183e-2 

-1.483700 

_sin(9.37°)_ 

yields 


y 

1.487 

sinOL 

:-2416_ 

(6-21) 


(6-22) 


which  implies  that  a  =  13.98°  Furthermore,  it  follows  that 


(z-Tw)^  •  (6-23) 

where  r„  =  3.3166  cm  (dome  radius  of  curvature).  This  defines  a 
circular  pupil  on  the  dome  of  radius  r  =  1.9  cm  (z  =  0.595  cm), 

computed  by  the  following  approximate  formula; 

r  .  tga.(t^-  z)  ,  .  (6.24) 
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The  quantities  a,  t^,  z  and  R^p  are  defined  in  Fig.  6-4. 


Figure  6-4.  Ray  tracing  variable  nomenclature. 


C.  RESULTS  AND  DISCUSSION 

The  characterization  of  fixed  pattern  noise  is  an  important 
aspect  in  evaluating  the  performance  of  a  staring  array.  Fixed 
pattern  noise  is  a  generic  term  used  to  describe  the  variations  in 
gain  and  offset  levels  observed  over  the  FPA  detector  elements. 

The  motivation  of  our  study  is  the  characterization  of 
fixedpattern  noise  induced  by  dome  emission. 

Computations  were  performed  for  the  case  of  a  hemispherical 
dome  corresponding  to  the  case  I  nose  configuration  (see  page  49)  of 
a  missile  at  Mach  4  at  sea  level  with  zero  incidence  angle.  Two 
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different  computations  were  performed,  the  first  considering  that 
the  dome  temperature  field  was  that  obtained  by  solving  the  thin- 
layer  Navier-Stokes  equations  and  considering  an  adiabatic  wall  as 
one  of  the  boundary  conditions;  in  the  other  case  the  dome 
temperature  was  considered  to  be  the  mean  value  of  the  first  case. 


DOME  TEMPERATURE 


Figure  6-5 .Thin-Layer-Navier-Stokes  computed 
dome  temperature. 

The  optical  system  was  considered  to  be  looking  straight  ahead 
and  just  the  portion  of  the  window  corresponding  to  the  system's 
window  pupil  was  considered.  The  collecting  lens  was  designed  using 
the  ZEMAX  XE  ray  tracing  program.  The  lens  data  are  given  in 
Appendix  E.  The  dome  material  chosen  was  CVD  Diamond.  The  waveband 
considered  was  from  S.O-lO.Sym.  Fig.  6-5  shows  the  dome  temperature 


128 


VBiTsus  off  axis  angls.  For  ths  cass  consid©r©d  th©  t©itvp©irat'ur©  field 
has  cylindrical  symmetry.  Fig.  6-6  shows  the  photon  irradiance  on  a 
64x64  FPA (Focal  Plane  Array)  due  to  this  temperature  distribution. 
Fig.  6-7  shows  the  photon  irradiance  on  a  64x64  FPA  based  on  the 
mean  dome  temperature  of  the  same  temperature  distribution.  The  mean 
dome  temperature  considered  is  1024  K. 

For  the  computations  cell  separation  of  SOpm  and  a  fill  factor 
of  1.0  are  considered.  A  fill  factor  of  1.0  although  unrealistic 
serve  the  purpose  of  showing  the  variation  of  the  photons  flux 
through  the  array. 

The  usual  means  of  characterizing  fixed-pattern  noise  is  by 
specifying  the  variance,  a,  of  the  pixel  levels, 

52  E  [l(x.ry^)  -  M]^ 

i -  .  (6-25) 

number  of  pixels 

The  variance  computation  is  often  written (schematically)  as: 

.  ~(l~Mp  -  P  -  (6-26) 

the  array  data  being  denoted  by  I(x,y)  and  the  mean  by  M.  A 
relatively  larger  variance  is  seen  for  arrays  exhibiting  an 
appreciable  amount  of  fixed-pattern  noise.  However,  a  specification 
of  variance  alone  gives  no  indication  as  to  how  the  non-uniformities 
are  distributed  over  the  array,  essentially  assuming  that  there  is 
an  equal  contribution  from  the  fixed-pattern  noise  at  all  spatial 
frequencies . 

Figure  6-6  shows  the  percentage  variation  of  the  quantum 
irradiance  on  a  64x64  focal  plane  array  due  to  a  dome  with  a  non- 
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uniform  temperature  field  distribution.  The  values  were  computed  in 
relation  to  the  maximum  quantiim  irradiance  at  the  center  of  the 
array.  The  maximum  quantum  irradiance  is  9.924x10^^  photons/sec  per 
pixel  at  the  center  of  the  array.  The  minimum  quantum  irradiance  is 
9.839x10’-^  photons/sec  per  pixel  at  the  corners  of  the  array.  The 
computed  mean  quantum  irradiance  is  6.2829  photons/sec  per  pixel. 
For  this  case  the  maximima  computed  irradiance  variation  is  -0.8%  at 
the  corners  of  the  array. 

Figure  6-7  shows  the  percentage  variation  of  the  quantum 
irradiance  on  the  focal  plane  array  due  to  a  dome  with  a  uniform 
temperature  field  distribution.  The  values  were  computed  in  relation 
to  the  maximum  quantum  irradiance  at  the  center  of  the  array.  The 
minimiam  quantum  irradiance  is  xlO^^  photons /sec  per  pixel  at  the 
corners  of  the  array.  The  computed  mean  quantum  irradiance  is 
6.2829  photons /sec  per  pixel.  For  this  case  the  maximum  computed 
irradiance  variation  is  -0.8%  at  the  corners  of  the  array. 

This  is  five  times  more  than  the  variance  found  for  a  uniform 
temperature  dome  at  the  average  computed  temperature.  (Fig. 6-7). 
The  computed  variance  is  1 . 26xl0^^photons^  per  pixel  for  the  case  of 
a  dome  with  uniform  temperature. 

Fig.  6-8  shows  a  detectivity  contour  plot  for  the  case  of 
uniform  dome  temperature.  For  that  purpose  the  following  equation 
[29]  was  used. 


K 


mean 


he 


(6-27) 


where  is  the  background  photon  flux  per  pixel  element  area.  The 
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quantuun  efficiency,  rj,  is  0.7. 
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64x64  FPA  Flioton  Irradiance  %variaton 


Figure  6-6,.  Photon  irradiance  percentual  variation  for  a 
64x64  FPA  considering  a  non-uniform  dome  temperature. 


64x64  FPA  Rioton  Irradiance  %variaton 


Figure  6-7.  Photon  irradiance  percentual  variation  on  a 
64x64  FPA  considering  a  uniform  dome  temperature. 


Figure  6-8.  Detectivity  contour  on  a  64x64  FPA  induced 
by  the  fixed  pattern  of  the  dome  background  irradiance 
(non-uniform  dome  temperature) 


VII.  CONCLUSIONS 


As  expected,  flight  at  high  supersonic  speeds  has  a  dramatic 
effect  on  target  detection  range  degradation  of  IR  guided  missiles 
due  to  the  aerokinetic  heating  of  its  IR  window  dome.  The 
temperature  increase  of  IR  dome  increases  the  background  quantum 
flux  density  that  irradiates  the  IR  sensor  leading  to  an  decrease  of 
the  signal-to-noise  ratio  and  the  correspondent  decrease  of  the 
NETD  of  the  detection  system.  An  equation  was  derived  that  shows 
this  effect.  For  Mach  niimber  speeds  of  2  and  4  the  detection  range 
degradation  for  the  cases  considered  reached  values  around  37%  and 
97%  respectively  for  array  non-uniformity  of  1%.  This  means  that  at 
high  supersonic  Mach  numbers  the  seeker  HR  ("Imaging  Infra-Red") 
system  becomes  "blind"  due  to  dome  emission.  For  lower  array  non¬ 
uniformities,  i.e.  0.001  percent  the  detection  range  degradation  was 
smaller  in  the  range  12%  and  22%  for  Mach  2  and  4  respectively. 

Since  the  integration  times  of  the  array  are  much  smaller  than 
the  missile  guidance  system  time  constant  (over  100  times)  the 
photon  irradiance  due  to  dome  emission  on  the  FPA  originates  a 
fixed  pattern  noise  on  the  detector  responsible  for  a  detectivity 
non-uniformity  of  the  array  pixels.  As  shown  in  Chapter  5,  the  IR 
seeker  detection  range  has  a  strong  dependance  of  the  detection 
array  non-uniformity. 
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APPENDIX  A.  CALCULATION  OF  DOME  THICKNESS 


Any  viable  missile  dome  structure  must  be  able  to  withstand  the 
mechanical  stresses  induced  by  aerodynamic  pressure  as  well  as 
thermal  stresses  induced  by  aerodynamic  heating.  In  a  thermally  thin 
regime,  the  TSR  (Temperature  Shock  Resistance)  performance  will  be 
maximized  by  making  the  dome  as  thin  as  possible.  The  TSR  indicates 
the  range  of  Mach  number  variation  that  an  IR  dome  can  withstand 
without  a  good  chance  of  catastrophic  failure [10] .  This  temperature 
corresponds  to  the  maximum  difference  in  temperature  between  the 
stagnation  temperature  and  the  temperature  of  the  inner  wall  of  the 
dome  that  the  dome  can  withstand.  This  leads  to  an  "optimum"  dome 
thickness  equal  to  the  minimum  dome  thickness  required  to 
withstand  the  aerodynamic  loads .  Since  the  emittance  of  the  dome 
increases  with  its  thickness  the  minimum  dome  thickness  will 
minimize  dome  self-emission  effects  as  well.  An  accurate  prediction 
of  the  dome  emission  effects  necessitates  an  accurate  computation  of 
the  dome  temperature  and  of  its  emittance,  and  by  consequence  of  the 
dome  thickness. 


1.  Pressure  Induced  Stress 

In  a  supersonic  flight  environment,  altitude  and  velocity 
determine  the  pressure  load.  We  should  follow  the  procedure  outlined 
in  section  II. C  in  order  to  calculate  the  minimum  dome  thickness  to 
withstand  the  aerodynamic  load. 


Ap  - 


I 


(A-1) 
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where 
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Y  -  1 

2y 

-1 

Y  -  1 

p. 

i  2  j 

1  Y  .  ij 

■ 

Y  .  1 J 

(A-2) 


Consider  a  medium  altitude  flight  profile.  From  U.S.  Standard 
Atmosphere  Tables  one  has  for  an  altitude  of  3  km: 

Y=1.4 

p„=10.2  psi. 

Evaluating  Eq.(A-2)  gives  a  pressure  jump  along  the  stagnation 
streamline  of  21.07.  Substituting  values  in  Eq. (A-1) ,  gives  Ap  = 
204.7  psi. 

Let  us  consider  that  a  "floating  attachment"  scheme,  like  a 
polyamide  adhesive,  is  used  to  decouple  the  dome  from  the  missile 
body.  The  maximum  tensile  stresses  failure  criterion  leads  to  the 
following  expression 


oJSF 

— -  -  0.581 

Ap 


(A-3) 


for  the  minimum  dome  thickness,  ,  where  Of  refers  to  the  nominal 
flexure  strength  of  the  material,  and  SF  is  an  appropriate  safety 
factor  that  yields  the  minimum  thickness  required  to  prevent 
catastrophic  failure  of  the  dome.  Since  currently  available  IR 
transmitting  dome  materials  are  brittle  ceramics,  they  may  exhibit 
a  wide  range  of  fracture  probabilities  under  apparently  identical 
loadings,  which  suggest  safety  factors  of  at  least  four  (4)  if  is 
the  characteristic  flexural  strength  derived  from  a  Weibull 
statistical  analysis .  The  dome  thickness  required  for  withstanding 
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aerodynamic  loads  Ap  with  good  probability  of  survival  can  be 
obtained  from  the  following  equation: 


min 

R 


I.  75 


V 


(A-4) 


For  a  truncated  hemispherical  dome  with  R  =  3.3166  cm, 
CVD  diamond  as  dome  material, 

L  .  .  1. 75  (3.32xl0-^m) 


1.41xlO‘‘Pa 
1 . 5x10^ Pa 


.  5. 4x10 -^m 


and  using 


(A-5) 


2 .  Thermally  Induced  Stress 

As  shown  in  Chapter  1  the  maximum  tensile  stress  should  be 
given  by  an  expression  such  as  [10] 


a 


max 


-^OLE'AT  , 
TSF 


(A-6) 


where  a.  -  thermal  expansion  coefficient, 

TSF  =  Thermal  Stress  Factor, 

E'  =  biaxaial  elastic  modulus, 

AT  =  temperature  differential  across  the  dome  wall. 

The  stress  induced  by  the  sudden  aerodynamic  heating  of  the  dome 
depends  on  the  ability  of  the  dome  to  conduct  the  heat  flow  and 
ultimately  upon  the  heat-flow  pattern,  which  is  controlled  by  the 
Biot  number. 
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Bi  . 


hL 

~T 


(A-7) 


where  h  =  heat  transfer  coefficient  at  the  outer  surface, 

L  =  dome  thickness, 

k  =  thermal  conductivity  of  dome  material. 

The  "optimum"  dome  thickness,  L  =  is  determined  by 

aerostructural  requirements  and  should  be  as  thin  as  possible  to 
enhance  the  thermal  shock  resistance  and  reduce  self-emisssion.  In 
assessing  the  "ultimate"  TSR  of  an  IR  dome,  the  Biot  nuinber  should 
be  expressed  in  the  following  way 
where  h'^^  refers  to  a  1-cm  radius  dome. 

At  Mach=4 ,  ■  altitude=3km 


Si'  - 


0.908^ 

~kaf^ 


(A-8) 


From  Fig  A-1, 


(0.908)  (3.3166cm)^^^ 
(20Pfcm-^K-^)  (150)^^^ 


1x10* 


(A-9) 


where  R  is  in  centimeters,  k  in  watts  per  centimeter  per  kelvin,  and 
Of  is  in  megapascals.  Considering  the  hemispherical  dome  window  case, 
R  =  3.3166  cm  and  CVD  Diamond  as  window  material,  we  conclude  that 
the  response  of  the  diamond  to  transient  heating  will  be  in  the 
thermally  thin  mode  (Bi  <1) .  Indeed  this  material  exhibits  extremely 
low  Biot  numbers (Bi<0 . 001 )  that  minimize  the  temperature  gradients 
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across  the  thickness. 


The  thermal  shock  temperature  (TST)  of  the  minimum  thick 
temperature  missile  domes  with  Bi*  that  represents  the  Biot  number 
for  the  most  severe  environment  based  on  projected  tactical 
requirements . 

For  thermally  thin  domes  the  peak  temperature  differential 
between  the  outer  and  inner  wall  is  given  by 


('At;  “  (T  .  T.  )xBi  ,  Bi<l 

P  St  z  w 


(A-10) 


In  the  thermally  thin  regime  the  temperature  of  the  back  face 
may  get  hot  long  before  the  front  reaches  the  recove3ry  temperature. 


Heat  Transfer  Coefficient  at  Stagnation 


I 


Figure  A-1.  Heat  transfer  Coefficient  at  Stagnation. 
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APPENDIX  B.  PROGRAM  TO  EXTRACT  FLOW  VARIABLES 

The  following  program  written  in  FORTIIAN  was  used  to  extract 
from  the  g.save  file  (PL0T3D  format)  the  values  of  the  flow 
variables  of  Table  4-2.  It  uses  the  grid. in  file  (PL0T3D  format)  and 
the  g.save  file  as  input  files.  By  a  simple  alteration  of  the  source 
code  it  is  possible  to  read  from  the  g.save  computed  file  flow  field 
values  for  a  specific  range  of  I(streamwise  direction), 
J (circumferencial  direction)  and  K(radial  direction). 
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C  Extract  x^y^z  vs.  Temperature^  Normalized  Temperature, 

C  Coefficient  of  Pressure  from  grid  and  Q  files  in  PLOT3D  format, 
C 

C  XYPLOT  and  MATLAB  formats. 

C 

PARAMETER  (MI=300 , MJ=1 00 , MK=50 ) 

DIMENSION  XYZ  (MI,MJ,MK,  3)  ,Q(MI,MJ,MK,  5) 

DIMENSION  X(MI) ,Y(MI) ,Z(MI) ,T(MI) 

CHARACTER* 80  FILE 
C 

PARAMETER  (GAMMA=1 . 4 ) 

C 

C 

C  Read  grid  file. 

C 

WRITE (*,1) 

1  FORMAT ('  Enter  PLOT3D  grid  filename: 

READ(*,2)  FILE 

2  FORMAT (A) 

OPEN (UNIT=1,  FILE-FILE, STATUS-' OLD'  , FORM-'  FORMATTED'  ) 

READd,*)  NI,NJ,NK 

IF  (NI.GT.MI  .OR.  NJ.GT.MJ  .OR.  NK.GT.MK)  THEN 
WRITE  (*,  3) 

3  FORMAT('  Grid  is  too  big,  exceeds  internal  dimensions.') 
GOTO  100 

END  IF 

READd,  *)  (  (  (  (XYZ  (I,  J,  K,NX)  ,  1  =  1,  NI)  ,  J=1,NJ)  ,K=1,NK)  ,  NX=1,  3) 

CLOSE (UNIT-1) 

C 

C  Read  Q  file. 

C 

WRITE {*, 4) 

4  FORMAT ('  Enter  PLOT3D  Q  filename:  ' , $) 

READ(*,2)  FILE 

OPEN (UNIT-1 , FILE-FILE, STATUS-' OLD' , FORM-' FORMATTED' ) 

READd,*)  NID,NJD,NKD 

IF  (NID.NE.NI  .OR.  N JD . NE . N J  .OR.  NKD.NE.NK)  THEN 
WRITE (*, 5) 

5  FORMAT ('  Q  file  dimensions  don' ' t  match  grid.') 

GOTO  100 

ENDIF 

READ  (1,  *)  FSMACH,  ALPHA,  RE,  TIME 

READ (1, *)  (  ( ( (Q(I,  J, K,NX)  , 1  =  1, NI)  , J=1,NJ)  , K-1,  NK)  , NX-1, 5) 

CLOSE (UNIT-1) 

C 

C  Get  range  of  (i,j,k). 

C 

C  WRITE(*,11) 

C  11  FORMAT  ('  Enter  IS,  IE,  JS,JE,  KS,KE;  ',$) 

C  READ(*,*)  IS, IE, JS, JE, KS, KE 

C 
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WRITE(*,  12) 

12  FORMAT ('  Enter  XYPLOT  Wall  Temp  filename:  ',$) 

READ (*,2)  FILE 

^  Ot’EN(UNIT=l,FILE=FILE,  STATUS='UNKNOWN'  ,  FORM=' FORMATTED'  ) 

GMl  =  GAMMA-1 . 

RHOINF  =  1. 

CINF  =  1 . 

PINF  =  1 . /GAMMA 
V2INF  =  0.5* (FSMACH/CINF) **2 
QINF  =  RH0INF*V2INF 
C  DO  20  K  =  KS,KE 

C  DO  20  J  =  JS,JE 

1  =  1 
J  =  1  , 

C 

DO  20  K  =  1, 35 

RHO  =  Q(I, J,K, 1) 

U  =  Q(I,  J,K,2)  /RHO 

V  =  Q(I, J,K, 3) /RHO 

W  =  Q(I, J,K, 4) /RHO 

EO  =  Q{I,  J,K,  5)  /RHO 

^2  =  0.5*(U**2  +  V**2  +  W**2) 

El  =  EO  -  V2 

P  =  GMl*RHO*EI 

CP  =  (P  -  PINF) /QINF 

T(I)  =  (P/PINF) / (RHO/RHOINF) *520*5/9 . 

Tratio  =  (P/PINF) / (RHO/RHOINF) 

X(I)  =  XYZ(I,J,K,1) 

Y(I)  =  XYZ(I,J,K,2) 

Z(I)  =  XYZ{I,J,K,3) 

WRITE (1,*)  X(I) ,Y(I) ,Z(I),T(I),Tratio,Cp 
20  CONTINUE 
CLOSE (UNIT=1) 

C 

C 

100  CONTINUE 
STOP 
END 
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APPENDIX  C.  OVERFLOW  INPUT/OUTPUT  FILE  SPECIFICATIONS 


The  following  list  describes  the  input/output  files  and  the  data 
they  contain: 


Filename 

Unit 

I/O 

Description 

1  ISSESIllllllli 

1 

grid  file  (PLOT3D  format) 

fort.2 

2 

input  from  PEGASUS 

3 

Q  restart  file  (PLOT3D  format) 

H 

Q  save  file  (PLOT3D  format) 

q.bomb 

4 

0  bomb  file  (PLOT3D  format) 

resid.out 

7 

BBSS! 

flow  solver  residual  history 

fomo.out 

8 

force  and  moment  coefficient  history 

rpmin.out 

9 

minimum  density/pressure/gamma  history 

turb.out 

11 

turbulence  model  residual  history 

12 

species  continuity  residual  history 

10 

both 

PEGASUS  interpolation  and  boundary  point  indices 

10 

both 

intermediate  storage  of  run  parameters 

flowvar.n 

10 

both 

intermediate  storage  of  flow  variables 

10 

both 

intermediate  storage  of  grid  and  metrices 

Table  C-1.  OVERFLOW  Input/Output  files. 
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APPENDIX  D.  LOWTRAN  ATMOSPHERIC  PROPAGATION  MODEL 


The  Atmospheric  Transmittance/Radiance  Computer  Code  LOWTRAN  6 
calculates  both  the  atmospheric  transmittance  and  radiance 
throughout  the  infrared  spectrum ( 0 . 25  to  28.5  gm)  [28]  .  It  uses  a 
single  parameter  band  model  for  molecular  absorption,  and  includes 
the  effects  of  continuum  absorption,  molecular  scattering,  and 
aerosol  absorption.  Refraction  and  earth  curvature  are  included  in 
the  calculation  for  slant  paths [28]. 

Five  seasonal  models  and  the  1962  U.S.  Standard  Atmosphere  are 
provided  as  inputs  to  the  LOWTRAN  program  using  typical  altitude, 
pressure,  temperature,  water  vapor  density  and  ozone  density 
profiles  for  each.  The  five  other  atmospheric  profiles  consist  of 
(1)  Tropical  (15°  N)  ,  (2)  Midlatitude  Summer  (45°  N-July)  ,  (3) 
Midlatitude  Winter  (45°  N-January) ,  (4)  Subarctic  Summer  (60°  N)  ,  and 
(5)  Subarctic  Winter  (60°  N- January) .  Additionaly  the  user  may  input 
radiosonde  data  if  desired  to  make  calculations  for  a  specific 
profile . 

The  LOWTRAN  code  includes  three  boundary  layer  aerosol  types: 
Rural,  Urban,  and  Maritime.  Due  to  the  nature  of  the  scenario,  only 
the  latter  was  used  to  simulate  an  engagement  of  an  anti-ship 
missile  in  open  waters.  These  aerosols  are  largely  sea-salt 
particles  in  the  lower  boundary  layer  caused  by  evaporation  of  sea- 
spray  droplets.  Together  with  a  background  aerosol  of  more  or  less 
pronounced  continental  characteristics  they  form  a  fairly  uniform 
maritime  aerosol  which  is  representative  of  the  lower  2  to  3 
kilometers  of  the  atmosphere  over  the  oceans . 


149 


These  aerosols  should  not  be  confused  with  the  heavy  direct 
sea-spray  aerosols  found  in  the  lower  10  to  20  meters  above  the 
ocean  surface  which  are  strongly  dependent  on  recent  wind  and  sea 
surface  turbulence.  The  extinction  and  absorption  coefficients  are 
calculated  as  a  function  of  the  relative  humidity  and  based  heavily 
upon  the  refractive  index  for  the  given  altitude  and  weather 
conditions . 

The  ability  to  calculate  slant-path  transmittance  is 
incorporated  in  the  LOWTRAN  VI  utilizing  a  data  base  of  seasonal  and 
latitude  dependencies  of  vertical  aerosol  distributions.  This  allows 
transmittance  calculations  of  medium  and  high  altitude  missile 
attack  profiles. 

The  LOWTRAN  code  calculates  transmittance  as  a  function  of 
"absorber  density"  for  the  path,  the  pressure,  and  temperature,  for 
the  particular  wavelength  band  chosen.  It  utilizes  both  empirical 
laboratory  data  and  available  molecular  line  constant  data  in 
performing  its  calculations . 

The  aerosol  population  found  over  the  world's  ocean  is 
significantly  different  in  composition  and  distribution  from  that  of 
a  continental  origin.  These  aerosols  are  largely  derived  from  the 
sea.  They  are  produced  by  the  evaporation  of  sea  spray  and  from  jet 
and  film  droplets.  Jet  droplets  are  ejected  into  the  air  by  the 
bursting  of  small  air  bubbles  at  the  sea  surface.  The  bursting  of 
the  bubble  film  leaves  behind  many  smaller  film  droplets  that  may 
also  be  diffused  into  the  air.  These  mechanisms  are  wind  dependent 
and  require  whitewater  phenomenon  to  produce  aerosol. 

Once  the  aerosol  droplets  are  airborne,  they  undergo  additional 
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sorting  and  mixing  processes .  The  marine  boundary  layer  is  usually 
capped  by  a  temperature  inversion  and,  within  the  boundary  layer, 
the  smaller  marine  aerosol  together  with  any  background  aerosol  form 
a  fairly  uniform  aerosol  spatial  distribution.  Once  introduced  into 
the  atmosphere  the  lifetime  of  an  aerosol  particle  is  dependent  on 
the  size  of  the  particular  aerosol  particle.  Those  with  very  small 
sizes  have  a  very  long  residence  time  in  the  boundary  layer  if  there 
are  no  washout  processes  taking  place.  On  the  other  hand,  those  with 
very  large  sizes  have  a  short  residence  time  and  do  not  contribute 
to  the  stationary  long-term  aerosol  population.  The  reader  is 
recommended  to  Reference [28]  for  more  details. 


APPENDIX  E.  SINGLE  LENS  DESIGN 


The  single  lens  used  on  the  radiative  transfer  of  the 
computation  of  detector  irradiance  was  designed  using  the  program 
ZEMAX-XE.  In  this  appendix  are  included  the  output  plots  of  such  a 
computation . 


PRRRXIRL  FQCRL  CHROMRTXC  SHIFT 


SEEKER  COLLECTING  LENS 
NON  NOV  J6  1992 

MRXIMUN  FOCRL  SHIFT  RRNGE :  91.5WS6  MICRONS 

DIFFRRCTION  LIMITED  RRNCE :  378.0376  MICRONS 


Figure  E-L  Seeker  collecting  lens  focal  chromatic  shift. 


Figure  E-2.  Seeker  collecting  lens  geometric  encircled  energy. 


Figure  E-3.  Seeker  collecting  lens  paraxial  chromatic  shift 
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APPENDIX  F.  "POLE-FIT"  EMISSIVITY  MODEL 


The  "classical  pole-fit"  model"  represents  the  complex 
permittivity,  e,  in  the  form  of  the  lattice  vibrations [24 ] : 


e(v ,T)  .  e 


Ae.  (T)Vj  (T) 


Vj (T)  - 


i (V,T)V 


(F-1) 


where  Ae^  ,  Yj  are  the  j'^'^  mode  strength,  line  width,  and  long  wave¬ 
length  transverse  optical  frequency,  respectively.  This  model 
adequately  describe  the  reflectance,  R,  in  the  infrared  region. 
Using  Eq.(F-l),  the  static  dielectric  constant,  e^iT)  =  e(0,T  ), 
becomes 


e^(T)  -  ejT)  .  ^Ae.(T)  . 

j 


(F-2) 


The  frequency  dependence  of  Yj  represents  a  cutoff  beyond  the 
one-phonon  region  caused  by  anharmonicities  of  the  potential.  The 
cutoff  is  empirically  represented  as  follows 


y.  .  y.  (T)  .  1  / 


.exp 


■ 

(  \ 

■ 

-a. 

-  I 

1 

V 


(F-3) 


where  oc  and  v  are  arbritary  parameters .  Michael  Thomas  team  at  APL 
has  been  most  successful  in  this  effort  by  choosing  a  =  4  and  ^  as 
1.1  times  the  highest  infrared  allowed  transverse  optical  optical 
mode  frequency. 
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As  Eq. (F-1)  suggests,  the  reflectance  is  a  function  not  only  of 
frequency  but  also  of  temperature.  The  temperature  dependence  of 
the  transverse  mode  frequencies,  ,  and  the  mode  strengths,  Ag^, 
can  be  represented  by  linear  equations  of  the  form 

(T)  -  (TJ  -  a.  (T  -  TJ  (P.4) 


and 


Ae^  (T)  .  Ae^  (TJ  .  b^(T  -  TJ 


(F-5) 


where  is  a  reference  temperature  and  a.^  and  are  mode-and 
material-specific  constants  coefficients.  The  dependence  of  the 
linewidth  on  temperature  can  be  described  by  a  quadratic  equation  in 
T  as 


Y-  Y- 

— (T)  ■  — ^  (T)  .  CAT  -  T  )  .  d.(T  -  T 

V  V  .  °  ^  o  J  O 

TOj  TOj 


(F-6) 


where  is  again  a  reference  temperature  and  and  dj  are  mode 
and  material-specific  constant  coefficients.  The  high-frequency 
peimittivity ,  ,  is  also  temperature  dependent  and  is  represented 
by 


G_  (T)  .  (TJ  .  e^(T  .  TJ  .  (p.7) 

ej=de^/ T=2n^  (dn„/dT) ;  thus,  in  most  cases,  e^-  is  determined  from 
dn/dT  measurements  at  the  helium-neon  laser  frequency. 
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APEENDIX  G.  RESULTS  FROM  ELECTROMAGNETIC  THEORY 

For  optically  smooth  surfaces  that  are  physically  and  chemically 
clean,  electromagnetic  theory  with  the  aid  of  Kirchoff's  law 
provides  a  prediction  of  the  magnitude  and  distribution  of  the 
directional  emissivity [12] .  The  radiation  emission  properties  of  the 
aforementioned  ideal  surfaces  are  sometimes  distinguished  by  using 
the  term  emissivity  instead  of  emittance. 

Radiant  energy  either  incident  on  or  leaving  from  a  surface 
can  be  resolved  into  two  components  of  polarization.  One  of  these 
components  is  polarized  parallel  to  the  plane  of  incidence  and  the 
other  is  polarized  perpendicular  to  the  plane  of  incidence.  The 
electromagnetic  theory  predicts  the  specular  reflectance  ff 
corresponding  to  each  component  of  polarization,  defined  as  the 
ratio  of  the  monochromatic  intensities  of  that  component  in  the 
reflected  and  incident  beams. 

The  expressions  for  p/  (perpendicular  polarized  component)  and 
P||®  (parallel  polarized  component)  constitute  the  well-known  Fresnel 
law  of  reflection  assuming  that  the  medium  adjacent  to  the 
reflecting  surface  has  a  refractive  index  of  unity  and  an  extinction 
coefficient  of  zero  (vacuum  and  gases) .  These  expressions  are 
derived  for  monochromatic  radiation.  However,  the  subscript  A  is 
omitted  in  order  to  avoid  complicating  an  already  cumbersome 
notation.  The  resulting  equations  are 

s  -  2a  cos  6  .  cos^  9 
a^  -  .  2a  cos  0  .  cos^  0 
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(G-lb) 


p,*re;  . 


-  2a  sin  6  tan  9  .  sin^  9  tan^  9 
a  ^  •  2a  sin  9  tan  9  .  sin^  9  taa^  9 


The  angle  0  is  the  specular-ray  direction  corresponding  to  an 
incident  beam  whose  angle  of  inclination  is  !P  =  6  .  The  quantities 
a  and  b  are  related  to  the  angle  of  reflection  0  and  to  the 
refractive  index  n  and  the  extinction  coefficient  k  of  the  material 
on  which  the  radiation  is  incident. 

2a^  .  ^  (n^  -  -  sin^  d)^  .  4n^k^  •  (n^  -  -  sin^  Q)  (G-2a) 

2b^  -  yj  (n^  .  .  sin^  9;^  .  4n^K^  -  (n^  -  k\.  sin^  d)  (G-2b) 


Figure  G-1.  Predictions  of  the  monochromatic  emissivity  and  the 
specular  reflectance  from  electromagnetic  theory. 

The  complex  index  of  refraction,  n,  is  defined  taking  the 
optical  constants  n  and  k  together,  as  follows 

n  .  n  -  ik  (G-3) 
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If  the  incident  radiation  is  unpolarized,  then  the  two  components 
of  polarization  in  the  incident  beam  are  of  equal  intensity.  In  this 
case  it  is  readily  verified  that  the  monochromatic  specular 
reflectance  is  the  average  of  p/  and  pf ,  that  is, 


(G-4) 


The  monochromatic  directional  absorptivity,  a{Q) ,  can  be  determined 
using  Kirchoff's  law. 


ol(Q)  .  1  -  pUQ)  ■  1  -  -/■p/('©;  •  p,^re;; 


(G-5) 


Due  to  the  fact  that  Kirchoff's  law  applies  directionally  and 
monochromatically  for  each  component  of  polarization,  one  can  write 
the  components 


zjQ)  -  ajQ)  .  1  -  p^re; 


(G-6b) 


2,  (Q)  .  a,  re;  -  I  -  pfre;  (G-6a) 

Furthermore,  the  directional  emissivity  of  the  mixed  radiation  is 
the  average  of  Sj_  and  ,  that  is, 

£  re;  -  ^  /■£,  re;  .  s,  re; ;  (g-?) 

Although  the  subscript  X  is  ommited  for  notational  simplicity,  the 
foregoing  equations,  as  well  as  those  that  follow,  are  to  be 
regarded  as  applying  monochromatically. 

The  predictions  of  electromagnetic  theory  for  Pj_®  and  pj|®  can  be 
employed  to  provide  an  expression  for  the  directional  emissivity. 


159 


Upon  substituting  equations  (G-1)  into  equations  (G-6)  e  (G-7) ,  one 
obtains 


see;  ■  -ejQ)  i 

2  ^ 


*  sin^Q 

cos^Q  (a^  ♦  *  2a  slnB  tanB  .  sin^Q  tan^Q 


(G8) 


where 


- 


4a  cosQ 

.  2a  sinQ  .  cos^Q 


(G-9) 


The  quantities  a  and  b  are  related  to  the  optical  constants  and  to 
the  angle  of  emission  6  by  equations  (G-2).  If  n  and  k  are  known, 
Eq. (G-8)  provides  a  prediction  of  the  directional  distribution  of 
the  monochromatic  emissivity  of  an  optically  smooth  surface. 

The  hemispherical  emissivity  e  can  be  determined  from  the 
directional  emittance  by  integrating  over  the  hemisphere.  By 
introducing  the  definitions  of  £■  and  the  emissivity  can  be 
computed  using  the  following  relation 

£  -  ^  £  (d)  cos  0  do)  . 

*ejn 

Eq. (G-8),  can  be  then  be  employed  to  provide  a  prediction  for  e.  In 
the  case  of  an  electric  nonconductor  (k  =  0) ,  the  execution  of  the 
integration  yields 
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,  ,  ,  (  \ 

'  2  3  '  3n  2(n^  ,  1)^  ""  i  n  - 

(n  .1)  (n^  .2n  -1)  4n^(n^  .1) 

.  1)^  (n  -1)  (n^  .  1)^  (n  -  1)^ 


in  which  is  the  emissivity  in  the  direction  normal  to  the  surface 


s 

n 


(G-12) 


The  ratio  of  the  hemispherical  to  the  normal  emissivity  ratio 
has  practical  significance  inasmuch  as  it  is  the  normal  emissivity 
that  is  more  often  determined  by  experiment,  whereas  it  is  the 
hemispherical  emissivity  that  is  usually  needed  in  the  computation 
of  the  radiant  interchange.  The  ratio  s/e^  is  presented  in  Fig.  G-2 
as  a  function  of  the  refractive  index  for  parametric  values  of  the 
extinction  coefficient  =  k/n.  Within  the  range  of  n  and  k  values 
that  correspond  to  real  materials,  the  £/e„  ratio  is  greater  than 
one  for  metals  and  less  than  one  for  dielectrics .  For  the  latter 
materials,  s  is  only  slightly  less  then 
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I  10  100 


REFRACTIVE  INDEX,  n 

Figure  G-2.  Ratio  of  the  hemispherical  to  the  normal 
emissivity  as  predicted  by  electromagnetic  theory. 
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APEENDIX  H.  RADIATION  TRANSFER  SOURCE  CODE 


The  program  listed  below  was  written  in  MATLAB  and  computes  the 
background  photon  flux  on  an  nxn  FPA  originated  by  a  hemispherical 
dome.  The  method  used  follows  the  procedure  outlined  in  the  section 
A  and  B  of  this  chapter.  The  program  assumes  that  the  dome 
temperature  has  cylindrical  symmetry  and  that  the  optical  system  is 
looking  straight  ahead. 

1 .  Input  Parameters 

The  temperature  distribution,  in  kelvins,  is  entered  in  the 
vector  tempdata.  The  first  element  is  the  dome  stagnation 
temperature . The  other  temperature  information  is  chosen  in  points 
downstream  in  order  to  ensure  that  a  reasonable  coverage  of  the  dome 
is  achieved.  The  coord  has  the  radial  coordinates  (measured  from  the 
dome  axis  in  cm)  of  the  points  corresponding  to  the  data  of  vector 
tempdata . 

Following  is  a  list  of  the  input  parameters: 

->  NW  -  radial  discretization  number  of  elements  of  the 
dome ; 

->  MW  -  number  of  sectors  used  to  discretize  the  dome; 

->  Rpupil  -  radius  of  the  pupil  measured  from  the  dome 
axis (m) ; 

->  R1  -  radius  of  curvature  of  lens  front  surface (m); 

->  R2  -  radius  of  curvature  of  lens  back  surf ace (m) ; 

->  Rdome  -  radius  of  curvature  of  the  dome(m); 

->  t  -  lens  thickness; 

->  tl  -  lens  separation  from  the  dome; 
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->  t3  -  backsurface  focal  length; 

->  lambdamin  -  cut-on  wavelength (microns ) ; 

->  lambdamax  -  cut-off  wavelength (microns ) ; 

->  ND  -  sqrt (number  of  array  pixels); 

->  cell  separation (microns ) ; 

->  em  -  dome  emissivity; 

->  n  -  index  of  refraction  of  the  media  suurronding  the 
lens  ; 

->  nl  -  index  of  refraction  of  the  lens; 

->  toptics  -  cold  optics  optical  transmission. 
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%%%%%%%%%%%.%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  @  LT  Clelio  D.  Ferreira  Leite  % 
%  Monterey,  16  Nov  1992  % 
%  % 
%  This  program  written  in  MATLAB  % 
%  computes  the  photon  irradiance  on  a  FPA  % 
%  originated  by  a  hot  window  hemispherical  dome  % 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear 

format  long  e 
ttiine=clock 

%%%%%%%%%%  input  parameters  %%%%%%%%% 

coord=[0  0.0078  0.0115  0.0132  0.0158  0.0175  0.0189]; 

Tempdata=[1222  1183  1150  1112  975  945  905]; 

Tempcoef =polyf it ( coord , Tempdata , 6 ) ; 

NW=41;%odd  # 

MW=60;%even  # 

Rpupil=,016; 

Rls=.600; 

R2=-.500; 

RW=. 03166 ; 
t=.003; 
tl=.01; 
t3=. 090999; 
lambdamin=8 . 0 ; 
lambdamax=10 . 5 ; 

ND=63;%even 
cellsep=30e-6; 
em= . 05 ; 
n=1.0; 
nl=4. 00431; 
toptics= . 5 ; 

%%%  Compute  coordinates  on  conjugate  window  and  detector  site  %%%%%%%%%%%%% 

xdlim=ND*cell sep/2 ; 

lambdamean=( lambdamax+lambdamin ) /2 ; 

Adbin=cellsep^2 ; 

Tstag=polyval ( Tempcoef , 0 ) %1144 ; 

2l=n*t/nl; 

22=l/R2*(nl/n-l)*(n*t/nl) ; 

2  3=1; 

zc=-zl/( ( z2+23 ) ) ; 

for  i=l:NW; 

r ( i ) =Rpupi 1* i/NW ; 

2(i)=sqrt(RW^2-r(i) ^2)+tl ; 

end 

for  i=l:2:NW-2; 

T ( i ) =polyval ( Tempcoef , r ( i+1 ) ) ; 

end 

for  i=l;NW; 

zl=z(i)+t, *z(i)/Rl*(n/nl-l)+n*t/nl; 

2  2=1/R2  *  ( nl/n- 1)  *  ( 2  ( i ) +t .  *  2  ( i ) /Rl*  ( n/nl-^l ) +n*t/nl ) ; 
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dywd=y de t ( n , m ) ; 
dzwd=t3-zc; 
normwdx=dxwd/rwd ; 
nonnwdy=dywd/rwd ; 
normwdz=d2wd/rwd ; 

dotprodw=inorm*nonawdx+ j  nonn*noniiwdy+knonn*noniiwd2 ; 
nonasurf=sqrt(  inonn^2+jnorm^2+knorin^2) ; 
norTnwd=sqrt  ( nonttwdx^2+nonnwdy^  2+nonnwdz  ^  2 ) ; 
costhetaw=abs ( dotprodw/ ( ( nonasurf  *nonnwd ) ) ) ; 
dotprodd=-normwd2 ; 
costhetad=abs ( dotpr odd/ ( normwd ) ) ; 

Qlainbda=ein*c3/(  lanibdamean^4*  ( exp(c2/(  lambdaniean*Tstag)  -1 ) ) )  *  C  laxabdamax-lambdainin 
Q=Qlaiiibda  *daw*Adbin*costhetad*cos  thetaw/rwd^  2  ; 

Qbin  ( n ,  m )  =Qbin  ( n  ,ia) +Q ; 
end 
end 
end 
end 

%%%  Compute  Photon  Flux  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
for  n=l:(ND+l)/2; 

elapsed_tiiiie=round( etime ( clock ,  ttime ) ) 
for  m=l:(ND+l)/2; 
for  k=l:2:NW-2; 
for  1=2:2:MW; 

% Compute  coordinates  on  conjugate  window  and  detector  site 
daw=pi*(r(k+2) . ^2-r(k) . ^2)/MW; 

rwd=sqrt(  (xconj(k+l,l)-xdet(n,m) )  .^2+(yconj(k+i,l)-ydet(n,m) )  ♦  ^2+(t3-zcon j(k+l,  1 

dxl=xcon j ( k+2 , 1 ) -xcon  j(k,l); 
dyl=yconj(k+2,l)-yconj(k,l) ; 
dzl=zconj(k+2,l)-2Conj(k,l) ; 

if  1=MW; 

dx2=xconj(k+l,l)-xconj(k+l ,MW-1) ; 
dy2=ycon j ( k+1 , 1 ) -ycon j ( k+1 , MW-1 ) ; 
d22=zconj(k+l,l )-zconj(k+l,HW-l) ; 
else 

dx2=xcon j (k+l, 1+1) -xcon j(k+l, 1-1) ; 
dy2=ycon j ( k+1 , 1+1 ) -ycon j (k+1 , 1-1 ) ; 
d22=2Conj(k+l,l+l)-2conj(k+l,l-l) ; 

end 

icross=dyl*dz2-d2l*dy2 ; 
jcross=-dxl*dz2+dzl*dx2 ; 
kcross=dxl*dy2-dyl*dx2 ; 

lcross=sqrt ( icross^  2+ jcross  ^  2+kcross  ^  2 ) ; 

inprm^icross/lcross ; 
jbp^=  jcross/lcross ; 
knorip«ckcross/lcross ; 

^^<J=xdet(n,in)-xconj(k+l  ,1) ; 

4y^d=ydet ( n , m ) -ycon j ( k+1 , 1 ) ; 
givd=t3-zcon j ( k+1 , 1 ) ; 

IWd=sqrt(dxwd^2+dywd^2+dzwd^2) ; 
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z3==nl/n*  ( z  ( i ) /Rl* ( n/nl-1 )  +n/nl ) ; 

2con(i)=-2l/( (z2+z3) ) ; 

Al=l+t/Rl* (n/nl-1 ) ; 

A2=1/R1* ( n/nl-l ) * ( 1+t/Rl * ( n/nl-1 ) +nl/n ) ; 

A=Al+zcon(i) .*A2; 
rconj ( i)=A*r( i) ; 

for  j=l:MW; 

xcon j ( i , j ) =rcon j ( i ) *cos ( 2*pi/MW*  j ) ; 
ycon j ( i , j ) =rcon j ( i ) *sin ( 2*pi/MW*  j ) ; 
zcon j ( i , j ) =abs ( 2Con( i ) ) ; 

end 

end 

x=xcon j; 
y=ycon j ; 
z=zconj; 

r=sqrt(x."^2+y,^2) ; 

xdet=zeros(ND) ; 
ydet=zeros ( ND ) ; 
for  n=l:ND+l 

for  in=l:ND+l; 

xdet  ( n  ,in)  =“Xdlim+ ( n-l )  *cellsep ; 
ydet  ( n ,  m )  =~xdl  im+  C  in-l )  *ce  llsep ; 

end 

end 

Qbin-zeros ( ND+1 ) 

axis ( 'square' ) 

Qlainbda=0 ; 

C2=l. 438864; 

C3=l. 8837627; 

%%%  Compute  photon  flux  originated  by  the  dome  central  element  %%%%%%%%%% 

%  rwd  -  virtual  distance  between  the  conjugate  surface  and  the  pixel  at 
%  detector ; 

%  daw  -  area  of  finite  element  at  the  conjugate  surface; 

%  costhetaw  -  angle  between  the  normal  to  the  conjugate  surface 
%  and  the  line  of  radiative  interchange; 

%  costhetad  -  angle  between  the  normal  to  the  detector 
%  and  the  line  of  radiative  interchange; 

%  Qlambda  -  photon  exittance  of  the  dome 

%  Q  -  photon  irradiance  originated  by  a  certain  portion  of  the  dome  in  a 
%  pixel  element; 

for  n=l:(ND+l)/2; 
for  m=l: (ND+l)/2; 

daw=pi *r ( 1 ) . ^  2 ; 

rwd=sqrt( (-xdet(n,m) ) . ^2+(-ydet(n,m) ) , ^2+(t3~zc) ^2) ; 

inorm=0 ; 
jnorm=0 ; 
knorm=l ; 

dxwd=xdet ( n , m ) ; 
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normwcix=dxwd/lwd ; 
nonnwdy=dywd/lwd ; 
nooiwdz=dzwd/lwd ; 

dotprodw=inorm*nonnwdx+  jnoritt*normwdy+)cnortn*norinwdz  ; 
nonasurf=sqrt  ( inorm^2+ jnonn^2+knorm^2 ) ; 
nonnwd=sqrt  ( nonawdx^  2+normwdy ^  2 +nonnwdz  ^  2 )  ; 
costhetaw=abs  ( dotprodw/ ( (normsurf  *norinwd) ) )  ; 
dotprodd=-nonnwdz ; 
costhetad=abs(dotprodd/(nonnwd) ) ; 

Qlainbda=em*c3/(  lambdamean^4* ( exp(c2/( lambdainean*T(k)  )-!)))*(  laxabdamax-lambdamin) 
Q=Qlainbda*daw*Adbin*costhetad*costhetaw/rwd^  2 ; 

Qbin(n,m)=Qbin(n,m)+Q; 

end 

end 

end 

end 

for  n=l:(ND+l)/2; 

for  itiia=(KD+l)/2+l:ND+l; 

Qbin(n,iniQ)=Qbin(n,ND+2-itmi) ; 
end 

end 

for  ia=l:ND+l; 

for  nn=(ND+l)/2+l2ND+l; 

Qbin ( nn , m ) =Qbin ( ND+2-nn , ra ) ; 
end 

end 

save  winl  Qbin 
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